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Abstract: 

In statistical QCD, color deconfinement and the properties of the quark-gluon plasma 
determine the in-medium behavior of heavy quark bound states. In high energy nuclear 
collisions, charmonia probe the partonic medium produced in the early stages of the 
interaction. We survey the present theoretical status and provide a critical evaluation 
of the charmonium production measurements in experiments at the CERN-SPS and the 
BNL-RHIC. 
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Introduction 



Statistical QCD predicts that strongly interacting matter undergoes a deconfining transi- 
tion to a new state, the quark-gluon plasma, in which the colored quarks and gluons are 
no longer bound to colorless hadrons. How can this state be investigated - which phenom- 
ena provide information about its thermal properties? The main probes considered so far 
are electromagnetic signals (real or virtual photons), heavy flavor mesons {Qq/Qq states), 
quarkonia {QQ pairs), and jets (energetic partons). In theory, the ultimate aim is to carry 
out ab initio calculations of the in-medium behavior of these probes in finite temperature 
QCD. In high energy nuclear collisions, we want to study the deconfinement transition 
and the QGP in the laboratory, to compare experimental results to QCD predictions. In 
our report, we want to consider as specific case the spectral analysis of quarkonia in a hot 
QGP and its application to nuclear collisions. 

The theoretical basis for this analysis is: 

• The QGP consists of deconfined color charges, so that the binding of a QQ pair in 
such a medium is subject to the effects of color screening. 

• The color screening radius ro(T) decreases with temperature T. 

• When td(T) falls below the binding radius r« of a QQ state i, the Q and the Q can 
no longer bind, so that quarkonium i cannot exist. 

• The quarkonium dissociation points Tj, specified through ro(Tj) ~ r i; determine 
energy density and temperature of the QGP, as schematically illustrated in Fig. 1. 
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Figure 1: Charmonium survival probabilities vs. energy density 

Clearly the essential aspect here is a precise determination of the dissociation points of 
the different quarkonium states, in temperature and/or energy density. This will therefore 
be the primary task of theoretical studies of the in-medium behavior of quarkonia. 

Quarkonium production in high energy collisions provides an excellent experimental tool 
to address a variety of challenging questions. 

• In pA collisions, quarkonium production can be studied in cold nuclear matter, 
in order to provide a reference to determine which features of AA collisions are 
"anomalous", i.e., due to a newly produced medium. 

• In AA collisions, quarkonium production can be measured as function of collision 
energy, centrality, transverse momentum, and A, providing a variety of different 
conditions of the produced medium. 
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• The onset of anomalous suppression for the different quarkonium states can be 
determined and correlated to thermodynamic variables, such as the temperature or 
the energy density. 

• The resulting thresholds in the survival probabilities Si of states % can then be 
compared to the relevant QCD predictions of Fig. I. 

In this way we could, at least in principle, obtain a direct comparison between experi- 
mental results and quantitative predictions from finite temperature QCD. 

Starting its operation at CERN in 1986, immediately after the theoretical prediction [1] of 
J/ip suppression as deconfinement signal, the experiment NA38 pioneered the systematic 
study of J/ip production and suppression. It was followed, some years later, by experi- 
ments NA50 and NA51, with adapted and upgraded versions of the same detector. The 
first of these experiments, NA38, had submitted a Letter of Intent in 1984 aiming at the 
study of thermal muon pair production. It was on the floor in 1986, but it had changed 
its main physics goal, as stated in its original proposal, to the study of J/tp suppression. 
The three experiments made use of an exceptionally performant muon pair spectrometer 
with good measurement resolution and high luminosity operating capability; both features 
were inherited from the previous NA10 experiment, which was designed for and devoted 
to the study of high mass Drell-Yan muon pairs. It took almost 20 years of efforts to 
improve the detectors and beams, refine the analyses, overcome various difficulties and 
collect useful data under various experimental conditions. This allowed double, and some- 
times triple checks of the outcoming results. In the meantime, new results are presently 
being provided by much more recent experiments, PHENIX at the BNL-RHIC collider 
and NA60 at the CERN-SPS. 

The main issue to be addressed is whether an anomalous J/ip suppression has or has 
not been discovered, as claimed in past publications. In order to answer this question, 
independently from any theoretical model, a mandatory prerequisite is to unambiguously 
define and precisely establish what is or should be considered normal J/ip behavior. We 
will see that the investment required to fill this part of the paradigm is always huge and, 
moreover, specific to each set of experimental conditions. It is, nevertheless, a sine qua 
non condition to avoid an insurmountable potential dead end to the problem. 



I Theory 

1 Heavy Quarks and Quarkonia 

Heavy quarks first showed up in the discovery of the J/t/j meson [2], of mass of 3.1 GeV; 
it is a bound state of a charm quark (c) and its antiquark (c) , each having a mass of some 
1.2-1.5 GeV. On the next level there is the T meson [3], with a mass of about 9.5 GeV, 
made up of a bottom or quark- antiquark pair (bb), with each quark here having a mass 
around 4.5 - 4.8 GeV. Both charm and bottom quarks can of course also bind with normal 
light quarks, giving rise to open charm (D) and open beauty (B) mesons. The lightest of 
these 'light-heavy' mesons have masses of about 1.9 GeV and 5.3 GeV, respectively. 
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The bound states of a heavy quark Q and its antiquark Q are generally referred to as 
quarkonia. Besides the initially discovered vector ground states J/if> and T, both the 
cc and the bb systems give rise to a number of other stable bound states of different 
quantum numbers. They are stable in the sense that their mass is less than that of two 
light-heavy mesons, so that strong decays into open charm or beauty are forbidden. The 
measured stable charmonium spectrum contains the IS* scalar i] c and vector J /if}, three 
IP states Xc (scalar, vector and tensor), and the 2S vector state if)', whose mass is just 
below the open charm threshold. There are further charmonium states above the if)'; these 
can decay into DD pairs, and we shall here restrict our considerations only to quarkonia 
stable under strong interactions. These are summarized in tables 1 and 2 for charmonia 
and bottomonia, respectively. The binding energies AE listed there are the differences 
between the quarkonium masses and the open charm or beauty threshold, respectively. 



state 


Vc 


J/i> 


XcO 


Xcl 


Xc2 


if;' 


mass [GeV] 


2.98 


3.10 


3.42 


3.51 


3.56 


3.69 


AE [GeV] 


0.75 


0.64 


0.32 


0.22 


0.18 


0.05 



Table 1: Charmonium states and binding energies 



state 


T 


Xbo 


Xbi 


Xb2 




Xbo 


X'bi 


Xbi 




mass [GeV] 


9.46 


9.86 


9.89 


9.91 


10.02 


10.23 


10.26 


10.27 


10.36 


AE [GeV] 


1.10 


0.70 


0.67 


0.64 


0.53 


0.34 


0.30 


0.29 


0.20 



Table 2: Bottomonium states and binding energies 



Quarkonia are rather unusual hadrons. The masses of the light hadrons, in particular 
those of the non-strange mesons and baryons, arise almost entirely from the interaction 
energy of their nearly massless quark constituents. In contrast, the quarkonium masses 
are largely determined by the bare charm and bottom quark masses. These large quark 
masses allow a very straightforward calculation of many basic quarkonium properties, 
using non-relativistic potential theory. It is found that the ground states and the lower 
excitation levels of quarkonia are very much smaller than the normal hadrons, and that 
they are very tightly bound. Now deconfinement is a matter of scales: when the sepa- 
ration between normal hadrons becomes much less than the size of these hadrons, they 
melt to form the quark-gluon plasma. What happens at this point to the much smaller 
quarkonia? When do they become dissociated? As already indicated, the disappearance 
of quarkonia can signal the presence of a deconfined medium [1], and the melting of spe- 
cific quarkonium states can determine the temperature of this medium [4]. Thus the 
study of the quarkonium spectrum in a given medium is akin to the spectral analysis of 
stars [5], where the presence or absence of specific excitation or absorption lines allows a 
determination of the temperature of the stellar matter. 

After these introductory comments, we shall review the basic properties of quarkonia and 
the dynamics of their dissociation, and then survey the different approaches to quarko- 
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nium binding in QCD thermodynamics. Following this, we consider the main theoretical 
features of quarkonium production in elementary as well as p — A and nucleus-nucleus 
collisions. With the theoretical basis provided, we then address the present experimental 
results. 

2 Quarkonium Binding and Dissociation 

We had defined quarkonia as bound states of heavy quarks which are stable under strong 
decay into open heavy flavor, i.e., M c5 < 2Md for charmonia and M b i < 2Mb for bottomo- 
nia. Since the quarks are so heavy, quarkonium spectroscopy can be studied quite well in 
non-relativistic potential theory [7] . The simplest ( "Cornell" ) confining potential [8] for 
a QQ at separation distance r has the form 

V(r)=ar-- (1) 

r 

with a string tension er ~ 0.2 GeV 2 and a Coulomb-like term with a gauge coupling 
a ~ 7r/12. The corresponding Schrodinger equation 

f 2 m c _ —V 2 + V(r)\ $,(r) = M^(r) (2) 
I m c J 

then determines the bound state masses Mj and the wave functions $j(r), and with 

( r 2) = J d 3 r r 2 |^(r)|7 J d 3 r |^(r)| 2 . (3) 

the latter in turn provide the (squared) average bound state diameters. 

The solution of eq. (2) gives in fact a very good account of the full (spin-averaged) 
quarkonium spectroscopy, as seen in Table 3. The line labelled AM shows the differences 
between the experimental and the calculated values; they are in all cases less than 1 %. 
Here r gives the QQ separation for the state in question, and the input parameters are 
m c = 1.25 GeV, m b = 4.65 GeV, ^ = 0.445 GeV, a = tt/12 [9]. 



state 


J/1> 


Xc 




T 


Xb 




Xb 




mass [GeV] 


3.10 


3.53 


3.68 


9.46 


9.99 


10.02 


10.26 


10.36 


AE [GeV] 


0.64 


0.20 


0.05 


1.10 


0.67 


0.54 


0.31 


0.20 


AM [GeV] 


0.02 


-0.03 


0.03 


0.06 


-0.06 


-0.06 


-0.08 


-0.07 


r [fm] 


0.50 


0.72 


0.90 


0.28 


0.44 


0.56 


0.68 


0.78 



Table 3: Quarkonium Spectroscopy from Non-Relativistic Potential Theory [9] 

We thus see that in particular the J / ifj and the lower- lying bottomonium states are very 
tightly bound {2M D B — M 3> Aq CD ) and of very small spatial size (r Cr^ 1.5 — 2 fm); 
here r h denotes the typical hadron diameter. Through what kind of interaction dynamics 
can they then be dissociated? 
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As illustration, we consider the collision dissociation of a J/ip. It is very small (rj/$ ~ 0.5 
fm) and hence can only be resolved by a sufficiently hard probe [10]. It is moreover tightly 
bound {2Mb — Mj/^ ~ 0.6 GeV), so that only a sufficiently energetic projectile can break 
the binding. In the collision with a normal hadron, the J/ip can thus only be dissociated 
through the interaction with a hard gluonic contituent of the hadron, not with the hadron 
as a whole (see Fig. 2). 




Figure 2: J/ip- hadron interaction 

The gluon momentum distribution g(x) in a hadron is determined in deep inelastic lepton- 
hadron scattering; with kh denoting the gluon momentum, x = kh/ph specifies the fraction 
of the incident hadron momentum pt carried by the gluon. For mesons, one finds for large 
momenta 

g(x)~(l-x) 3 , (4) 
so that the average momentum of a hadronic gluon is 

<*>» = \(Ph)- (5) 
For thermal hadrons in confined matter, (p^) ~ 3T, with T < 175 MeV, so that with 

(k) h = -T < 0.1 GeV < AE ~ 0.6 GeV (6) 

5 

the gluon momentum is far too low to allow a dissociation of the J/ip. 

On the other hand, the average momentum of a deconfined thermal gluon in a quark-gluon 
plasma will be 

(k g ) ~3T (7) 

and for T > 1.2 T c , this provides enough energy to overcome the J/ip binding. We thus 
expect that a hot deconfined medium can lead to J/ip dissociation, while the gluons 
available in a confined medium are too soft to allow this. 

To make these considerations quantitative, one first has to calculate the cross-section for 
the gluon-dissociation of a J/ip, a QCD analogue of the photo-effect. This can be carried 
out using the operator product expansion [10,11], and the result is 

1 (fc/Ag, - l) 3 / 2 
~ ml (k/AE^ (8) 
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with AEj/,p = 2Md — Mj/$. The corresponding cross-section for the hadron dissociation 
is then obtained by convoluting the gluon-dissociation cross-section (8) with the hadronic 
gluon distribution function g(x), which for J/^-meson interactions leads to 

(Jh-JW ^ Ogcom(l - A /A) 5 - 5 (9) 

with A ~ (s - Mfy/M^ and A ~ (M h + AEy). Here ff gcom ~ vr(r J/ ^/2) 2 ~ 2 mb is 
the geometric J/ip cross-section and denotes the mass of the incident meson. In Fig. 
3, we compare the two dissociation cross-sections (8) and (9) as function of the incident 
projectile momentum. 




k [GeV] 

Figure 3: Gluon and hadron J/ip dissociation cross-sections [10] 

This result confirms the qualitative arguments given at the beginning of this section: 
typical thermal gluon momenta near 1 GeV produce a large dissociation cross-section, 
whereas hadron momenta in a thermal range (up to 2 - 3 GeV) still lead to a vanishingly 
small cross-section. In other words, the J/ ip should survive in confining media, but become 
dissociated in a hot quark-gluon plasma [10]. 

The calculations leading to eq. (8) are exact in the large quark mass limit niq — > oo. It 
is not clear if the charm quark mass really satisfies this condition, and so the dissociation 
cross section in J/^-hadron collisions has been discussed in other approaches, some of 
which lead to much larger values [12,13]. The basic question is apparently whether or not 
the charmonium wave function has enough overlap with that of the usual hadrons to lead 
to a sizeable effect. It may well be that this question is resolvable only experimentally, by 
testing if slow charmonia in normal nuclear matter suffer significant dissociation. Such 
experiments are definitely possible [14]. 

3 Thermal Quarkonium Dissociation 

We shall here address the first of our "basic" questions: Is it possible to specify the state 
of strongly interacting matter in terms of observable quantities calculated in QCD? 
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3.1 Interaction Range and Color Screening 



Consider a color-singlet bound state of a heavy quark Q and its antiquark Q, put into 
the medium in such a way that we can measure the energy of the system as function of 
the QQ separation r (see Fig. 4). The quarks are assumed to be heavy so that they are 
static and any energy changes indicate changes in the binding energy. We consider first 
the case of vanishing baryon density; at T = 0, the box is therefore empty. 



Q Q Q Q 

r 

4 ► 




Figure 4: String breaking for a QQ system 



In vacuum, i.e., at T = 0, the free energy of the QQ pair is assumed to have the string 
form [8] 

F(r) ~ or (10) 

where a ~ 0.16 GeV 2 is the string tension as determined in the spectroscopy of heavy 
quark resonances (charmonium and bottomonium states). Thus F(r) increases with sep- 
aration distance; but when it reaches the value of a pair of dressed light quarks (about 
the mass of a p meson), it becomes energetically favorable to produce a qq pair from the 
vacuum, break the string and form two light-heavy mesons (Qq) and (Qq). These can 
now be separated arbitrarily far without changing the energy of the system (Fig. 4). 

The string breaking energy for charm quarks is found to be 

F = 2(M D - m c ) ~ 1.2 GeV; (11) 
for bottom quarks, one obtains the same value, 

F = 2(M B - m b ) ~ 1.2 GeV, (12) 

using in both cases the quark mass values obtained in the solution leading to Table 3. 
Hence the onset of string breaking is evidently a property of the vacuum as a medium. It 
occurs when the two heavy quarks are separated by a distance 

r ~ 1.2 GeV/cr ~ 1.5 fm, (13) 

independent of the mass of the (heavy) quarks connected by the string. 

If we heat the system to get T > 0, the medium begins to contain light mesons, and the 
large distance QQ free energy F(oo,T) decreases, since we can use these light hadrons 
to achieve an earlier string breaking through a kind of flip-flop recoupling of quark con- 
stituents [15], resulting in an effective screening of the interquark force (see Fig. 5). Near 
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Figure 5: In- medium string breaking through recoupling 

the deconfmement point, the hadron density increases rapidly, and hence the recoupling 
dissociation becomes much more effective, causing a considerable decrease of F(oo,T). 

A further increase of T will eventually bring the medium to the deconfmement point 
T c , at which chiral symmetry restoration causes a rather abrupt drop of the light quark 
dressing (equivalently, of the constituent quark mass), increasing strongly the density of 
constituents. As a consequence, F(oo,T) now continues to drop sharply. Above T c , light 
quarks and gluons become deconfined color charges, and this quark-gluon plasma leads 
to a color screening, which limits the range of the strong interaction. The color screening 
radius ro, which determines this range, is inversely proportional to the density of charges, 
so that it decreases with increasing temperature. As a result, the QQ interaction becomes 
more and more short-ranged. 

In summary, starting from T = 0, the QQ probe first tests vacuum string breaking, then 
a screening-like dissociation through recoupling of constituent quarks, and finally genuine 
color screening. In Fig. 6, we show the behavior obtained in full two-flavor QCD for the 
color-singlet QQ free energy as a function of r for different T [16]. 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 



Figure 6: The color singlet QQ free energy F(r,T) vs. r at different T [16] 

It is evident in Fig. 6 that the asymptotic value F(oo,T), i.e., the energy needed to 
separate the QQ pair, decreases with increasing temperature, as does the separation 
distance at which "the string breaks". Here we consider the latter to be defined by 
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the point beyond which the free energy remains constant within errors. The behavior 
of both quantities is shown in Fig. 7. Deconfmement is thus reflected very clearly in 
the temperature behavior of the free energy of a heavy quark potential: both the string 
breaking energy and the interaction range drop sharply around T c . The latter decreases 
from hadronic size in the confinement region to much smaller values in the deconfined 
medium, where color screening is operative. 




Figure 7: Asymptotic free energy and interaction range at different temperatures 



The in-medium behavior of heavy quark bound states thus serves quite well as probe 
of the state of matter in QCD thermodynamics. We had so far just considered QQ 
bound states in general. Let us now turn to a specific state such as the J/ij). What 
happens when the range of the binding force becomes smaller than the radius of the 
state? Since the c and the c can now no longer see each other, the J/ip must dissociate for 
temperatures above this point. Hence the dissociation points of the different quarkonium 
states provide a way to measure the temperature of the medium. The effect had already 
been illustrated schematically in Fig. 1, showing how with increasing temperature the 
different charmonium states disappear sequentially as function of their binding strength; 
the most loosely bound state disappears first, the ground state last. 

Moreover, since finite temperature lattice QCD also provides the temperature dependence 
of the energy density, the melting of the different charmonia or bottomonia can be specified 
as well in terms of e. In Fig 8, we illustrate this, combining the the energy density 
calculated in lattice QCD [18, 19] and the force radii from Fig. 7. It is evident that 
although ip' and Xc are expected to melt around T c , the corresponding dissociation energy 
densities will presumably be quite different. 

To make these considerations quantitative, we thus have to find a way to determine the 
in-medium melting points of the different quarkonium states. This problem has been 
addressed in two different approaches: 

• Solve the Schrodinger equation (2) with a temperature-dependent potential V(r, T), 
obtained either through model considerations or from heavy quark lattice studies. 

• Calculate the quarkonium spectrum directly in finite temperature lattice QCD. 

Clearly the last is the only model-independent way, and it will in the long run provide 
the decisive determination. However, the direct lattice study of charmonium spectra 
has become possible only quite recently, and corresponding studies for bottomonia are 
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0.5 1.0 1.5 T/T c 

Figure 8: Charmonium dissociation vs. temperature and energy density 

still more difficult. Hence much of what is known so far is based on Schrodinger equation 
studies with different model inputs. To illustrate the model-dependence of the dissociation 
parameters, we first cite early work modelling the temperature dependence of V(r,T), 
then some more recent studies based on lattice results for F(r,T), and finally summarize 
the present state of direct lattice calculations of charmonia in finite temperature media. 

3.2 Potential Model Studies 

The first quantitative studies of finite temperature charmonium dissociation [4, 5] were 
based on screening in the form obtained in one-dimensional QED, the so-called Schwinger 
model. The confining part of the Cornell potential (1), V(r) ~ err, is the solution of the 
Laplace equation in one space dimension. In this case, Debye-screening leads to [21] 

V(r,T)~arl— — — J, (14) 

where n(T) denotes the screening mass (inverse Debye radius) for the medium at tem- 
perature T. This form reproduces at least qualitatively the convergence to a finite large 
distance value V(oo,T) = cr//x(T), and since fi(T) increases with T, it also gives the 
expected decrease of the potential with increasing temperature. Combining this with the 
usual Debye screening for the 1/r part of eq. (1) then leads to 

V(r, T) ~ art 1 ' 6 - ] - -e~» r = - (l - e^ r ) - -e~» r (15) 
[ fj,r J r fx 1 ' r 

for the screened Cornell potential. In [4], the screening mass was assumed to have the 
form /x(T) ~ 4 T, as obtained in first lattice estimates of screening in high temperature 
SU (N) gauge theory. Solving the Schrodinger equation with these inputs, one found that 
both the if)' and the Xc become dissociated essentially at T ~ T c , while the J /if) persisted 
up to about 1.2 T c . Note that as function of the energy density e ~ T 4 , this meant that 
the J/ip really survives up to much higher e. 

This approach has two basic shortcomings: 

• The Schwinger form (14) corresponds to the screening of ar in one space dimension; 
the correct result in three space dimensions is different [21]. 
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• The screening mass /i(T) is assumed in its high energy form; lattice studies show 
today that its behavior near T c is quite different [23]. 

While the overall behavior of this approach provides some first insight into the problem, 
quantitative aspects require a more careful treatment. 

From thermodynamics it is known that F can be decomposed into the internal energy U 
and the entropy S of the system, 

F = U-TS. (16) 

Since 

F = -T\nZ (17) 
in terms of the partition function Z, we have 

u-r-^f) as) 

and 

*--(£)■ 

Given F(r, T), eqns. (18) and (19) determine U(r, T) and S(r, T). It should be emphasized 
that the F(r,T) calculated in finite temperature lattice QCD specifies the difference in 
free energy between a system with and a system without a QQ pair, so that also U and 
S are such differences. Since U (r, T) thus is the expectation value of the the Hamiltonian 
with a heavy quark pair, minus that without such a pair, it is the relevant quantity to 
determine the binding potential. We shall see shortly, however, that other effects enter 
as well. 

When the first lattice results for the color averaged value F av (r,T) of the free energy 
became available, they were used to determine melting points for the different charmonium 
states [22]. It was assumed that in the thermodynamic relation (18) the entropy term 
—T(dF/dT) could be neglected, thus equating binding potential and free energy, 

V(r,T)~F av (r,T). (20) 

Using this potential in the Schrodinger equation (2) specifies the temperature dependence 
of the different charmonium masses. On the other hand, the large distance limit of V(r, T) 
determines the temperature variation of the open charm meson D, 

2Md(T) ~ 2m c + V(oo, T) (21) 

In fig. 9, we compare the resulting open and hidden charm masses as function of tem- 
perature. It is seen that the ip' mass falls below 2 M D around 0.2 T c , that of the Xc at 
about 0.8 T c ; hence these states disappear by strong decay at the quoted temperatures. 
Only the ground state J/ip survives up to T c and perhaps slightly above; the lattice data 
available at the time did not extend above T c , so further predictions were not possible. 

The main shortcoming of this approach is also quite evident. The neglect of the entropy 
term in the potential reduces V(r,T) and hence the binding. As a result, the D mass 
drops faster with temperature than that of the charmonium states, and it is this effect 
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Figure 9: Temperature dependence of open and hidden charm masses [22] 



which leads to the early charmonium dissociation. Moreover, in the lattice studies used 
here, only the color averaged free energy was calculated, which leads to a further reduction 
of the binding force. 

Today, finite temperature lattice QCD with two light quark flavors [16, 17] provides the 
color singlet free energy F(r, T) of a static QQ pair in a strongly interacting medium at 
temperature T. Here F(r, T) is again the difference between a system with and a system 
without a QQ pair; it was shown in Fig. 6 for several temperatures below and above the 
deconfinement point T c . For very small separation r <C T^ 1 , the QQ pair as a color neutral 
entity does not "see" the medium. Hence the results for different temperatures must 
coincide in the short distance limit and fall on the T = curve F(r,T = 0) = or — a/r; 
this is used to normalize the results for T > T c . Given F(r,T), we can now extract the 
corresponding color singlet behavior of the internal energy and the entropy; in Fig. 10, we 
show schematically the resulting behavior of the different thermodynamic potentials as 
function of r at some fixed T > T c and as function of T for r — > oo [24] . At very short QQ 
separation distance and finite T, for r <C T~ x , the QQ feels no medium and the medium 
does not see the pair; hence U(r,T) = F(r,T) = F(r,0) and S(r,T) = 0. The onset of 
the entropy starts around r ~ r£>(T), i.e, when the QQ separation distance reaches the 
value of the screening length of the medium. 




Figure 10: Thermodynamic potentials as function of r and T 



For a quantitative potential theory study, we now have to extract the in-medium heavy 
quark potential V(r,T) from the color singlet internal energy obtained in lattice QCD. It 
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appears that so far this task is not solved in an unambiguous way, in spite of a number 
of attempts. We therefore summarize these and their results, indicating the remaining 
problems. 

The internal energy U is the expectation of the Hamiltonian 7i, in our case the difference 
in interaction energy between a system with and one without the heavy quark pair at rest 
in the medium. Hence we have 

U(r,T) = — + V(r,T)), (22) 
m 

where the first term gives the relative kinetic energy and V(r, T) the interaction potential. 
The lattice results we want to use are obtained for a static QQ pair, so that there is no 
relative kinetic energy, and hence U (r, T) is our effective potential. However, at very large 
separation, the Q and the Q do not see each other any more and hence cannot interact. 
Nevertheless, they will still polarize the medium in their vicinity and thus make U(oo,T) 
different from zero. This polarization cloud can be interpreted as an effective heavy quark 
thermal mass m(T) > m c . Near T c , the gluonic correlation length becomes very large (or 
even diverges), thus increasing the size of the clouds. It is this effect which causes the 
dramatic increase in U(oo,T) as T approaches T c (see the right part of Fig. 10). As the 
Q and the Q get closer to each other, the polarization clouds begin to overlap and hence 
interact. The interaction enhances the binding potential between the two heavy quark 
components, and this cloud-cloud interaction is what leads U(r,T) near td to overshoot 
the T = form. The free energy does not contain this component and hence approaches 
the Cornell potential. For a further discussion, see [25]. 

The relevant Schrodinger equation thus becomes 

2m c - — + U(r,T)\ $,(r,T) = M t (T)^(r,T), (23) 
m c J 

and the resulting dissociation temperatures are listed in Table 4. As we shall see shortly, 
they agree quite well with the values presently quoted in lattice studies. However, the 
inclusion of the cloud-cloud interaction has been put to question, and so the Schrodinger 
equation has also been solved with potentials of an intermediate form V(r,T;x) = 
x U(r,T) + (1 — x) F(r,T), with < x < 1. The results tend towards those of Ta- 
ble 4 for x — > 1 and give lower dissociation temperatures as x — > [26-28] . 



state 




Xc(lP) 


^(2S) 


T(1S) 


XbO-P) 


T(2S) 


X b (2P) 


T(3S) 


T d /T c 


2.10 


1.16 


1.12 


> 4.0 


1.76 


1.60 


1.19 


1.17 



Table 4: Quarkonium Dissociation Temperatures [9] 

We thus have to conclude that at this time, a final answer is still not provided in the 
potential study approach. Besides the mentioned potential ambiguities, a further diffi- 
culty shared by all potential models is that the dissociation points are defined as those 
temperature values for which the diameter of a given state i diverges (see Fig. 11). 
As seen, this leads to "bound state" regions in which r, ^> T _1 and in fact surpasses 
the normal hadronic size r^. One may expect [29] that under such conditions, thermal 
dissociation will break up any bound state, so that the actual bound state survival ends 
when its radius reaches the size of the screening radius [30]. 
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Figure 11: T- dependence of the bound state diameter for the J/if> [9] 



3.3 Charmonium Correlators 

The direct spectral analysis of charmonia in finite temperature lattice has come within 
reach only in very recent years [31-36]. It is possible now to evaluate the correlation 
functions G»(t, T) for charmonium quantum number channels % in terms of the Euclidean 
time t and the temperature T. These correlation functions are directly related to the 
corresponding spectral function <7j(M, T), 



which describe the distribution in mass M at temperature T for the channel in question. 
Here the kernel 



provides the relation betweeen the imaginary time r and the cc energy u. The inversion 
of eq. (24) with the help of the maximum entropy method (MEM) provides the desired 
spectrum [37]. In Fig. 12, schematic results at different temperatures are shown for the 
J I if> and the Xc channels. It is seen that the spectrum for the ground state J / if) remains 
essentially unchanged even at 1.5 T c . At 3 T c , however, it has disappeared; the remaining 
spectrum is that of the cc continuum of J /if) quantum numbers at that temperature. 
In contrast, the Xc is already absent at 1.1 T c , with only the corresponding continuum 
present. 

These results are very promising: they indicate that in a foreseeable future, the disso- 
ciation parameters of quarkonia can be determined ab initio in lattice QCD. After first 
calculations performed in quenched QCD [31-35], i.e., without dynamical quark loops, 
there now are also results from two-flavor QCD [36], and these fully support the late 
dissociation of the J /if). The widths of the observed spectral signals are at present de- 
termined by the precision of the lattice calculations; to study the actual physical widths, 
much higher precision is needed. Moreover, one has so far only first signals at a few 
selected points; a temperature scan also requires higher performance computational facil- 
ities. Since the next generation of computers, in the multi-Teraflops range, is presently 
going into operation, the next years should bring the desired results. 




(24) 



K(u,r,T) 




(25) 



sinh(u;/2T) 
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So far, in view of the mentioned uncertainties in both approaches, the results from direct 
lattice studies and those from the potential model calculations of the previous section 
appear compatible. We note at this point, however, that numerous attempts to obtain 
information on charmonium survival by comparing the calculated correlator to reference 
correlators constructed from models of low or zero temperature spectra and relating the 
results to potential studies [38-40] have so far not really provided unambiguous informa- 
tion. Hence further work on the theory of in-medium charmonium behavior, in lattice 
studies as well as in analytic approaches [41-45], is clearly called for. 

4 Charmonium Production in Hadronic Collisions 

4.1 Elementary Collisions 

The hadroproduction of charmonia occurs in two stages. The first stage is the production 
of a cc pair; because of the large quark mass, this process is well described by perturbative 
QCD (Fig. 13). A parton from the projectile interacts with one from the target; the 
(non-perturbative) parton distributions within the hadrons are determined empirically in 
other reactions, e.g. by deep inelastic lepton-hadron scattering. The produced cc pair is 
in general in a color octet state. In the second stage, it neutralises its color and then 
eventually forms physical resonances, such as J/ip, Xc or if)'. Color neutralisation occurs 
by interaction with the surrounding color field; this and the subsequent resonance binding 
are presumably of non-perturbative nature. 

On a fundamental theoretical level, color neutralization is not yet fully understood. How- 
ever, the color evaporation model [46] provides a simple and experimentally well-supported 
phenomenological approach. In the evaporation process, the c and the c can either com- 
bine with light quarks to form open charm mesons (D and D) or bind with each other to 
form a charmonium state. The basic quantity in this description is the total sub-threshold 
charm cross section R c5 , obtained by integrating the perturbative cc production cross sec- 
tion a over the mass interval from 2m c to 2m/> At high energy, the dominant part of R c5 
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(a) (b) 




(c) 

Figure 13: Lowest order diagrams for cc production in hadronic collisions, through gluon 
fusion (a,b) and quark- ant iquark annihilation (c). 

comes from gluon fusion (Fig. 13a), so that we have 

r2mo r 

Rcc(s) ^ ds dxidx 2 gp(xi) g T (x 2 ) a(s) 5(s-xix 2 s), (26) 

J2m c J 

with gp(x) and gr(x) denoting the gluon densities [47], X\ and x 2 the fractional momenta 
of gluons from projectile and target, respectively; a is the gg — > cc cross section. In 
pion-nucleon collisions, there are also significant quark-antiquark contributions (Fig. 13c), 
which become dominant at low energies. The basic statement of the color evaporation 
model is that the production cross section of any charmonium state % is a fixed fraction 
of the subthreshold charm cross section, 

<n(s) = ft R C c(s), (27) 

where fi is an energy-independent constant to be determined empirically. It follows 
that the energy dependence of the production cross section for any charmonium state 
is predicted to be that of the perturbatively calculated sub-threshold charm cross section. 
As a further consequence, the production ratios of different charmonium states 

^4 = k = const. (28) 
W) fi 

must be energy-independent. Both these predictions have been compared in detail to 
charmonium and bottomonium hadroproduction data over a wide range of energies [48]; 
they are found to be well supported, both in the energy dependence of the cross sections 
and in the constancy of the relative species abundances. Let us consider in more detail 
what this tells us about the hadronization of charm quarks. 

We recall that the relative abundances of light hadrons produced in hadron-hadron and 
e + e~ interactions follow the statistical pattern governed by phase space weights [49,50]: 



20 



the relative production rates are those predicted by an ideal resonance gas at the confine- 
ment/deconfinement transition temperature T c ~ 175 MeV. For two hadron species % and 
j that implies at all (high) collision energies 



for the ratio of the production rates, with di for the degeneracy (spin, isospin) and rrii for 
the mass of species i. For strange particles, the rates (29) overpredict the experimental 
data; this can, however, be accommodated by a common strangeness suppression factor 
7 S ~ 0.5 — 0.7, applied as 7™ if the produced hadron contains n strange quarks [51]. A 
recent alternative explanation for strangeness abundances is based on the assumption of 
a color event horizon with corresponding Unruh radiation [52]. 

For the hadroproduction of charm, such a statistical description does not work, as seen 
in three typical instances [48,53]: 

• The total cc cross section increases with energy by about a factor ten between 
y^s = 20 and 40 GeV, while the light hadron multiplicity only grows by about 20%. 
Hence the ratios of hadrons with and without charm are not energy-independent. 

• From p — p data one finds for J/ip production a weight factor fj/^ ~ 2.5 x 10~ 2 . 
Since the subthreshold cc cross section is about half of the single D production 
cross section [53], this implies R(j/^)/d — 10 -2 ; the ideal resonance gas gives with 
R{j/ip)/D — 10~ 3 a prediction an order of magnitude smaller. Of the total charm 
production, more goes into the hidden charm sector than statistically allowed. 

• For the production ratio of ip' to J/ip, which have the same charm quark infrastruc- 
ture, one finds experimentally over a wide range of collision energies R^i /(jm) — 0.23. 
This energy-independent ip' to J / ip ratio can be accounted for in terms of the char- 
monium masses and wave functions; it disagrees strongly with the statistical pre- 
diction, which gives with Rp/(j/$) — 0.045 a very much smaller value. The same 
holds true for the other measured charmonium states [9]. 

Charm production in elementary collisions thus does not seem to be of statistical nature. 
It appears to be determined by parton dynamics at an early stage [54], rather than by 
the phase space size at the confinement temperature. 

Although the color evaporation model provides a viable phenomenological description of 
the hadroproduction of quarkonia, leading to correct quantitative predictions up to the 
highest energies under consideration, it cannot predict the fractions fi of the hidden charm 
cross sections, and it can even less describe the space-time evolution of color neutralization. 
For charmonium production in p-A and A-B collisions, the latter is crucial, however, and 
hence a more detailed description of color neutralization is needed. 

In the first step of the collision process, a cc pair is formed through a hard process, with 
very short formation time r c5 ; to reach subsequently J / ip quantum numbers, the cc must 
be in a color octet state. The color octet model [55] proposes that this cc then combines 
with a soft collinear gluon to form a color singlet (cc—g) state. After a relaxation time r 8 , 
this pre-resonant (cc—g) state turns into the physical cc mode by absorbing the accompa- 
nying gluon, with similar formation processes for Xc and tp' production. The color octet 




(29) 
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Figure 14: The evolution of J/ip production 



model encounters difficulties if the collinear gluons are treated perturbatively, illustrating 
once more that color neutralization seems to require non-perturbative elements [56,57]. 
However, it does provide a conceptual basis for the evolution of the formation process 
(see Fig. 14). The color neutralization time r§ of the pre- resonant state can be esti- 
mated [58]; it is essentially determined by the lowest momentum possible for confined 
gluons, r 8 — (2m c AqcD) 1 ^ 2 — 0.25 fm. The resulting scales in J/ip formation are il- 
lustrated in Fig. 15. The formation time for the actual physical ground state J/ip is 
presumably somewhat larger than r 8 ; although rj/^/2 ~ r 8 , the heavy c quarks do not 
move with the velocity of light. For the larger higher excited states, the formation times 
will then be correspondingly larger. 
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Figure 15: Scales of J/ip production 

There is one further important feature to be noted for J/ip hadroproduction. The J ftp's 
actually measured in hadron-hadron collisions have three distinct origins: about 60 % 
are directly produced IS charmonium states, while about 30 % come from the decay 
X C (1P) — > J/ip+ anything, and the remaining 10 % from ip'(2S) — > J/ip+ anything [59]. 
Such feed-down also occurs in T production [60]. In all cases, the decay widths of the 
involved higher excited states are extremely small (less than one MeV), so that their life- 
times are very long. The presence of any medium in nuclear collisions would therefore 
affect these excited states themselves, and not their decay products. 

4.2 p —A Collisions 

In p —A collisions, the presence of normal nuclear matter can affect charmonium produc- 
tion, so that such collisions provide a tool to probe the effect of confined matter. Nuclear 
effects can arise during the entire evolution of J / ip production, and several different phe- 
nomena have been studied in considerable detail. We note in particular: 

• The presence of other nucleons in the nucleus can modify the initial state parton dis- 
tribution functions, which enter in the perturbative cc production process illustrated 
in Fig. 13. 
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• Once it is produced, the cc pair can be dissociated in the pre-resonance as well as 
in the resonance stage, due interactions with nucleons during its passage through 
the target nucleus. 

In both cases, the crucial quantity is the momentum of the charmonium state as measured 
in the nuclear target rest frame. 

Since we eventually want to probe the effect which the 'secondary' medium produced 
by nucleus-nucleus collisions has on charmonium production, it is of course essential to 
account correctly for any effects of the nuclear medium initially present. Let us therefore 
first summarize the main features observed for charmonium production in p — A collision 
experiments [62]. 

• At fixed collision energy, quarkonium production rates per target nucleon decrease 
with increasing A. 

• The production rates decrease for increasing J/ip momentum as measured in the 
nuclear target rest frame. 

• The nuclear reduction at p — N mid-rapidity appears to become weaker with in- 
creasing collision energy. 

• For fixed collision energy, mass number A and J/tp rapidity, the reduction appears 
to increase with the centrality of the collision. 

• At sufficiently high momentum in the target rest frame, the different charmonium 
states appear to suffer the same amount of reduction, while at lower energy, the 
ip' is affected more than the J ftp. 

At present, there does not seem to exist a theoretical scenario able to account quantita- 
tively for all these observations. In fact, so far not even a common scaling variable for 
the different effects has been found. Shadowing would suggest scaling in the fractional 
target parton momentum x 2 , while absorption of the pre-resonance state would point to 
the fractional beam momentum xf- Neither is in good accord with the data. The problem 
has recently been addressed in the context of parton saturation and color glass conden- 
sate formation [61]. Here we shall concentrate on indicating some operational methods 
to specify nuclear effects in p — A collisions in a way that can be extended to A — B 
collisions; for further discussions, we refer to reviews [62,63]. The problems encountered 
in formulating a theoretical description of charmonium production on nuclear targets un- 
derline again the crucial importance of having p — A data in order to arrive at a correct 
interpretation of A — A results. 

We now return in some more detail to the two main aspects arising for charmonium 
production in pA collisions, the (initial state) nuclear modification of the parton distri- 
bution function and the (final state) absorption of the nascent charmonium state during 
its traversal of the nucleus. 

Assuming again gluon fusion as the dominant high energy process of cc formation, we 
have for J/ip production in pA collisions (see eqs. (26) and (27)) 

r2m D r 

Rj/ip( s ) - fj/i> / d § / dx 1 dx 2 g p (x u ii 2 ) g A (x 2 ,iJ, 2 ) a(s) S(s - xix 2 s), (30) 
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where p denotes the scale of the probe in the evaluation of the gluon distribution functions 
in either the proton (g p ) or the nucleus (g A ) . The nuclear modification of the distribution 
functions have been studied in detail in different approaches [64]. In Fig. 16 we show 
a typical result for the ratio R A (x) of the gluon distribution function of a nucleon in a 
heavy nucleus, relative to that in a single nucleon. It is shown as function of the Bjorken 
variable x, specifying the fraction of the nucleon momentum carried by the gluon. We 
note that there are essentially four different regimes. Near x — 1, Fermi motion inside 
the nucleus leads to an enhancement; following this is a suppression (the EMC effect) due 
to nucleon-nucleon interactions. Around x ~ 0.1, we then again have an enhancement 
( "anti-shadowing" ) , followed by a suppression at very small x ( "shadowing" ) . High energy 
collisions generally access regions x < 0.1, so that anti-shadowing and shadowing are the 
phenomena of particular interest. Depending on the collision energy and the momentum 
of the observed J/ip, the production rate can thus be either enhanced or reduced by the 
initial state nuclear modification of the gluon distribution function. In order to determine 
the effect of the final state (cold nuclear matter or, in AA collisions, a possible QGP) on 
charmonia, such initial state modifications must evidently be brought under control. 
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Figure 16: Nuclear modifications of the gluon distribution function 

The second aspect is the dissociation of the nascent J/ip in cold nuclear matter; the cc 
pair has thus been formed, with whatever nuclear modifications of the parton distribution 
functions applied for its formation rate, and now traverses the remaining part of the 
nucleus. To specify the dissociation, we note [66] that according to the Glauber formalism 
of nuclear scattering theory, a cc pair formed at point z in a target nucleus A has a survival 
probability 

S? = J d 2 bdzp A (b, z) exp {-(A - 1) £° dz' p A (b, z>) ^ iss } , (31) 

where the integration covers the path, at impact parameter b, remaining from z out 
of the nucleus, and where a l diss describes the overall "absorption" effect on the observed 
charmonium state i along the path. The result is then averaged over impact parameter 
and path lengths. The traversed medium of nucleus A is parametrized through a Woods- 
Saxon density distribution p A (z), and by comparing Si to data for different targets A, 
the effective dissociation cross section can be obtained for the J/ip and ip' absorption in 
nuclear matter. The effect of the charmonium passage through the nucleus will arise from 
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a superposition of the different stages; but if part of the passage is carried out as physical 
resonance, higher excited states should lead to larger absorption cross-sections than the 
much smaller ground state J/ip. 

The way charmonium production in pA (or dA) collisions enters the analysis is thus 
clear: with the initial state modifications of the parton distribution functions taken into 
account, the dissociation cross section in cold nuclear matter has to be determined. With 
this given as a baseline, one can then turn to AA collisions and look for the effects of a 
newly produced medium. 

It should be noted here that in some studies [66] the initial and final state effects have 
been parametrized jointly in the form of a common "dissociation cross-section", which 
includes both sources of modification. With the help of this cross section, normal nuclear 
matter effects were then calculated for nucleus-nucleus collisions. Evidently such a pro- 
cedure makes sense only if neither of the nuclear modifications depend sensitively on the 
momentum of produced charmonium state, or if the accessible momentum range is the 
same and very small both in p — A and in nuclear collisions. For the SPS data, the latter 
seems to be the case. 

4.3 Nuclear Collisions 

The basic assumption in the attempt to create deconfined matter through nuclear col- 
lisions is that the excited vacuum left after the passage of the colliding nuclei forms a 
thermal medium. This picture is schematically illustrated in Fig. 17. A charmonium 
state produced in such a collision will in its early stages first be subject to the possible 
effects of the nuclear medium, just as it is in p — A collisions, and then, after the nuclei 
have separated, encounter the newly produced medium. 




before collision nuclear medium produced medium 

Figure 17: Collision stages 



The Glauber formalism used above to calculate the survival probability of an evolving 
charmonium state in a p — A collision now has to be extended to A — B interactions [66]. 
The survival probability at impact parameter b now becomes 

Sf B (b) = J d 2 s dz dz 1 p A (s, z)p B {b - s, z') x 

~ ^ / a dZA PA ( S ' ZA) a ' diss ~( B ~ l ) B dzB P B ^ h ~ S ' Zb ) a diss \ > (32) 
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as extension of Eq. (31). Here Zq specifies the formation point of the ccg within nucleus 
A, Zq its position in B. With the dissociation cross sections cr l diss determined in p — A 
collisions, Eq. (32) specifies the 'normal' survival probability, i.e., that due to only the 
nuclear medium. 

To use quarkonia as probes for the produced medium, we now have to study how the 
behavior observed in A — A collisions differs from this predicted pattern. Several possible 
and quite different effects have been considered as consequences of the produced medium 
on charmonium production. 

• Suppression by comover collisions: A charmonium state produced in a primary 
nucleon-nucleon collision can be dissociated through interactions with the con- 
stituents of any medium subsequently formed in the collision. Such dissociation 
could in principle occur in a confined [70] as well as in a deconfined medium [10]. 

• Suppression by color screening: If the produced medium is a hot QGP, it will dissoci- 
ate by color screening the charmonium states produced in primary nucleon-nucleon 
collisions [1]. Since the different states have different dissociation temperatures, a 
sequential suppression pattern will specify the relevant thresholds [4,5,71]. 

• Enhancement by regeneration: In the hadronization stage of the QGP, charmonium 
formation could occur through the binding of a c with a c from different nucleon- 
nucleon collisions. If the total number of available cc pairs considerably exceeds their 
thermal abundance, such statistical regeneration could enhance hidden relative to 
open charm production, as compared to hadron-hadron collisions [72-75]. 

In addition, the partonic initial state of the colliding nuclei, which leads to the formation 
of the produced medium and to that of charmonium states, will change its nature for 
large A and sufficiently high yfs; eventually, parton percolation (saturation or color glass 
formation) can lead to a very different medium, with possible effects on production and 
binding of charmonia [76]. 

Is it possible for experiment to distinguish between these different scenarios? Before 
turning to the experimental situation, we want to discuss in some more detail the salient 
features of each approach. 

4.3.1 Suppression by Comover Collisions 

If the charmonium state moves in a random scattering pattern through the produced 
medium, its survival rate is approximately given by 

Si = exp{— (XjnTo ln[n/n/]} (33) 

with <7j denoting the dissociation cross section, n the initial density of the medium after 
a formation time To, and rif the 'freeze-out' density, at which the interactions stop. 

Since the cross section for J/ ip break-up through gluon collisions is large [10] and the gluon 
density high, there will be significant charmonium suppression in a deconfined medium, 
even if this is not thermalized. In an equilibrium QGP, this dissociation is presumably 
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accounted for by color screening, provided the effect of the medium on the width of the 
surviving states is also calculated. 

Charmonium dissociation by interaction with hadronic comovers has received consider- 
able attention in the past [70]. However, if one restricts the possible densities to values 
appropriate to hadronic matter (n^<0.5 fm -3 ) and the cross sections to those obtained in 
section 3, the effect of hadron dissociation is negligible. Even a cross section increase to 
the high energy limit still leads to less than 10% effects. More recent analyses [77] thus 
conclude that a hadronic medium will not result in significant suppression. 

In Fig. 18 we illustrate schematically the overall behavior expected for J/ip dissociation 
through comover collisions, assuming that beyond a deconfinement threshold, the co- 
mover density increases with energy density in a monotonic fashion, with little or no prior 
suppression in the hadronic regime. 
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Figure 18: J/ip suppression by comover 
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Figure 19: Sequential J/ip suppression 
by color screening 



4.3.2 Suppression by Color Screening 

The theoretical basis of this effect has been considered in detail in chapters 3 and 4; 
the color field between the heavy quarks becomes modified due to the presence of a 
medium of unbound color charges. The results obtained for this effect in statistical QCD 
are as such model-independent, once all calculational constraints are removed. What 
is speculative and model-dependent is its application to nuclear collisions: these do not 
necessarily produce the medium studied in thermal QCD, and the different evolution 
stages in nuclear collisions can introduce factors not present in the study of equilibrium 
thermodynamics, such as the oversaturation of cc pairs just mentioned. 

If the medium produced in high energy nuclear collisions is indeed the quark-gluon plasma 
of statistical QCD, and if charmonium production can be treated as a distinct process 
within such an environment, then the effect of color screening seems clear. The partition- 
ing of the cc pairs produced in nucleon-nucleon collisions into hidden and open charm is 
non-statistical, favouring the hidden charm sector because of dynamical binding effects. 
Color screening destroys these and hence strongly suppresses charmonium production 
rates relative to those observed in elementary interactions. 
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A crucial feature of J/ip suppression by deconfinement is its sequential nature [5,30,71, 
78]. In the feed-down production of J/ip, the produced medium affects the intermediate 
excited states, so that with increasing temperature or energy density, first the J/^'s 
originating from ip' decay and then those from Xc decay should be dissociated. Only 
considerably higher temperatures would be able to remove the directly produced J/^'s. 
Such a stepwise onset of suppression, with specified threshold temperatures, is perhaps the 
most characteristic feature predicted for charmonium as well as bottomonium production 
in nuclear collisions. It is illustrated schematically in Fig. 19, where we have defined the 
J/ip production probability to be unity if the production rate suffers only the calculated 
nuclear suppression. 

4.3.3 Enhancement through Regeneration 

In charmonium hadroproduction, J/^'s are formed because some of the cc pairs produced 
in a given collision form a corresponding bound state. In a collective medium formed 
through superposition of many nucleon-nucleon (NN) collisions, such as a quark-gluon 
plasma, a c from one NN collision can in principle also bind with a c from another 
NN collision. This 'exogamous' charmonium formation at a later evolution stage could 
lead to enhanced J/ip production, provided the overall charm density of the medium at 
hadronization is sufficiently high and provided the binding force between charm quarks 
from different sources is large enough [72-74]. 

The production of cc pairs in primary collisions is a hard process and thus grows in A —A 
interactions with the number of nucleon-nucleon collisions; in contrast, the multiplicity 
of light hadrons grows roughly as the number of participant nucleons. Hence the relative 
abundance of charm to non-charm quarks will be higher in A — A than in p —p collisions. 
Moreover, the cc production cross section increases faster with energy than that for light 
hadron production. The two effects together imply that in the medium produced in 
energetic A —A collisions, the ratio of charm to non-charm quarks is initially much higher 
than in a equilibrated QGP. Whether or not this results in a charmonium regeneration 
depends on two factors. On one hand, the initial charm oversaturation must be preserved, 
so that the total charm abundance in non-thermal. On the other hand, it is necessary 
that the binding potential of random pairs into charmonia is sufficiently strong. 

Charmonium regeneration in nuclear collisions will be addressed in a separate section 
of this Handbook [75]. We only note here that several studies have led to considerable 
enhancement factors [75,79], in the strongest form even predicting a large overall en- 
hancement of J/ip production in A — A collisions relative to p —p results scaled by binary 
collisions. One crucial prediction of the approach is the increase of the enhancement with 
centrality, as shown in Fig. 20, because of the corresponding increase in the number of 
collisions and hence of the number of cc pairs. Another is that the distributions of the 
observed charmonia in transverse momentum as well as in rapidity must be given as con- 
volutions of the open charm distributions [80]. What this means for transverse momentum 
spectra will be addressed in the next subsection. 
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Figure 20: J /if) enhancement by statistical regeneration 
4.4 Transverse Momentum Behavior 

The production pattern of charmonia as function of their transverse momentum can pro- 
vide information about the production process, the evolution of the resonance formation 
and of that of the produced medium. We begin with the production process. 

The transverse momentum distribution of charmonia measured in pA as well as in AA 
collisions is generally broadened in comparison to that in pp interactions. The main 
effect causing this is the collision broadening of the incident gluons which fuse to make 
cc pairs. A standard random walk analysis gives for the squared transverse momentum 
of the produced J/ip [81,82] 

(Pt)pA=(Pt) pp + N?6 (34) 

in p— A and to 

(p 2 t)aa = (P 2 T ) PP + N AA 5 (35) 

in A — A collisions. Here N A denotes the average number of collisions of a projectile par- 
ton in the target nucleus A, and N AA the sum of the average number of collisions of a 
projectile parton in the target and vice versa, at the given centrality. The parameter 5q 
specifies the average "kick" which the incident parton receives in each subsequent colli- 
sion. The crucial parameters are thus the elementary (pt)pp from p—p interactions and the 
value of 5q, determined by corresponding p— A data; both depend on the collision energy. 
The ^-dependence of N A as well as the behavior of N AA as function of centrality can be 
obtained through a Glauber analysis including absorption in cold nuclear matter, thus 
specifying the "normal" centrality dependence of (pt)aa- in the absence of any anoma- 
lous supppression, this would be the expected behavior of the average J/ip transverse 
momentum. 

Given the sequential suppression by color screening, the J/ip's observed for energy densi- 
ties in a range above the onset of anomalous suppression, e(2S), e(lP) < e < e(lS) are the 
directly produced IS states unaffected by the presence of the QGP. They should there- 
fore still show the normal broadening pattern (35), increasing linearly with the number of 
collisions as long as e < e(lS) [71]. This broadening is a memory of the initial state and 
hence essentially absent if the J/ip is formed only at the hadronization point, where such 
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memory has been destroyed. Thus J/ip production through regeneration should show a 
flat distribution as function of the number of collisions [80], in contrast to the rise ex- 
pected in the sequential suppression scenario. The behavior of (pt)aa thus should provide 
a clear indication of how the observed J/^'s were produced. In Fig. 21 we illustrate the 
different patterns expected. 
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Figure 21: Transverse momentum broadening in sequential suppression vs. statistical 
regeneration scenarios 

A further interesting signal is how an anomalous suppression behaves as function of px- 
If the nascent J / ip is formed as a small color singlet gradually expanding to its physical 
size, then it can be out of the deconfining medium either spatially [83] or temporally [84] 
before suffering serious suppression, provided it has a high enough transverse momentum. 
The "end" of anomalous suppression as function of pt could thus give indication on the 
size or life-time of the QGP. 

More recently, calculations have appeared [85] which are based on the AdS/CFT cor- 
respondence applied to QCD; they indicate the opposite effect, suggesting that J/^'s of 
high transverse momentum should suffer more suppression. It must be kept in mind, how- 
ever, that binding and the dissociation of specific bound states are clearly non-conformal 
phenomena, so that such modelling is not necessarily valid. 

5 Summary of the Theoretical Status 

The theoretical work of the past years has done much to clarify the questions that must 
be answered in order to reach final conclusions. 

In statistical QCD, improved MEM techniques, better statistics (from faster computers) 
as well as larger lattices (on larger computers), promise that in the foreseeable future we 
will have reliable first principle calculations of the in-medium behavior of quarkonia. 

These will hopefully then be related to the results from heavy quark studies in finite 
temperature lattice QCD. In this work (see e.g. [16,17,86]), there are at present still several 
unsolved questions, concerning in particular gauge invariance and the specification of 
specific color states, and the possible in-medium transitions between color states. Recent 
analytical work may provide help in dealing with them [42-45]. 
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In parallel, potential theory studies are expected to provide further insight into the melting 
of quarkonia in a hot QGP. One of the basic issues to be addressed here is if and how the 
multi-component problem near T c can be addressed. In particular, the role of the gluonic 
dressing in heavy quark binding has to be taken into account [87]. 

Once the in-medium behavior of quarkonia is clarified in statistical QCD, one can hope to 
address the behavior observed in high energy nuclear collisions. If these indeed produce 
a QGP in the sense of statistical QCD, sequential melting threshold could provide a 
direct quantitative connection between theory and experiment. If, on the other hand, 
the measured quarkonia are largely due to regeneration at hadronization, i.e., to the 
combination of secondary Q and Q from different nucleon-nucleon interactions, it is not 
evident how such a connection can be obtained. It would, however, provide clear evidence 
for the production of a deconfined thermal medium. 



II Experiment 

1 Charmonium Experiments at the CERN-SPS 

At the CERN-SPS, experiments NA38, NA50 and NA51 have systematically measured 
charmonium production with incident protons, Oxygen, Sulphur and Lead beams at var- 
ious incident momenta. The latest version of the detector [88] was based on a muon 
pair spectrometer able to stand high intensity incident beam fluxes. The spectrometer 
was further equipped with three devices leading to three independent estimates of the 
centrality of the collision. The Electromagnetic Calorimeter measured E^, the integrated 
flux of neutral transverse energy, mainly due to neutral pions produced in the collision. 
The Zero Degree Calorimeter measured Eznc, the forward energy carried by the beam 
spectator nucleons. The Multiplicity Detector counted the charged particles produced in 
the reaction. All these three measurements are related to the impact parameter of the 
collision. 

The properties of the detectors allowed the collection of large samples of data with incident 
protons and ions. This, in turn, made possible the study of both J / ip and ip' production 
with minimal systematic uncertainties, thanks to the simultaneous detection of Drell-Yan 
events which could be used as appropriate "debiasing" tools. 

From the very beginning, it was taken for granted that the study of J / ip production from 
different nuclear targets would provide, by some kind of extrapolation, the appropriate 
reference baseline, i.e., the "normal" charmonium behavior, relative to which the heavy 
ion collision specific features would be easily identifiable. Ideally, the reference data 
would have had to be collected with proton beams of the same energy of that of the 
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corresponding ion beams, namely 200 GeV, the energy of the Oxygen and Sulphur beams, 
and later 158 GeV, the energy of the Lead and Indium beams. For low cross-section 
measurements like charmonium production, it was highly desirable that proton beams be 
primary beams, directly extracted from the accelerator. Because of various reasons, low 
energy primary proton beams were not available at CERN. It took 18 years since the start 
of the heavy ion program until they finally became available, for about 3 days, and could 
thus be used by experiment NA60. 

1.1 The Nuclear Dependence of Charmonium Production 

The uncontroversial feature observed since long in charmonium production measurements 
in p — A reactions is that, at fixed collision energy, quarkonium production rates per 
target nucleon decrease with increasing A, the target atomic mass number. This effect 
has been traditionally, and somewhat empirically, quantified with the parametrization 
&pA = &o x A a which leads to measured values of a lower than unity. More recently, it has 
been assumed that some kind of cc state is created in the reaction which is dissociated, or 
rather prevented to finally form a bound state, through interaction with the surrounding 
nuclear matter. In other words, the cc pairs are created proportionally to the number 
of nucleon-nucleon collisions or, equivalently, to A. Some interact with the surrounding 
nuclear medium while moving through the target and are no more able to reach a bound 
state. The final number of charmonium bound states reaching the detector exhibits 
therefore an apparent suppression with respect to the original A dependence of the cc 
pair production. An appropriate variable to parametrize the measured J / ip cross-sections 
should then be the number of nucleons that the created state can potentially interact with 
[89]. This number can be calculated from the product p x L, where p is the nuclear density 
distribution and L the length of nuclear matter the cc state traverses while escaping 
from the interaction region. The charmonium survival probability can then be calculated 
as a function of an "absorption" or dissociation cross-section <7 a b s , using the Glauber 
formalism of nuclear scattering theory or, alternatively, a simplified " pL" exponential 
parametrization. Comparison with the experimentally measurable survival probability 
allows to determine <7 a b s • As a matter of fact, systematic experiments can even provide 
some guidance for the validity of this elaborated view of nuclear effects in charmonium 
production. 

Let us underline here that the a parametrization can be used to fit a set of several 
measurements made under same conditions with different nuclear targets. The procedure 
then becomes an hypothesis test of such a parametrization. It can also be assumed 
that nature behaves according to this law and apply the parametrization to only two 
different targets. The procedure then would just provide an estimate of the numerical 
value of a, taking for granted that such a parametrization correctly accounts for the 
elementary process. Identical procedures can be followed under the dissociation cross- 
section interpretation and resulting parametrization. Finally, let us point out that a and 
the more elaborated (x abs are obviously correlated. 
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1.2 Normal Charmonium Production 



The first significant samples of charmonium events produced in p — A reactions were 
collected by experiment NA38 in 1988, with 450 GeV incident protons on various targets. 
Later on, they were complemented with measurements obtained, as a by-product, from 
experiment NA51 which provided the data allowing to extract the charmonium production 
cross-sections for p-p and p-d collisions at 450 GeV [90] . All these data, reanalyzed, with 
identical procedures as those used to analyze the heavy ion data, led to the production 
cross-sections plotted in Fig. 1 as a function of the atomic mass number of the target [91]. 

Some time later, this set of measurements could be complemented with other results from 
the same experiment NA38, both with protons and with Oxygen and Sulphur beams in- 
cident on Copper and Uranium targets, all for collisions at 200 GeV per nucleon. When 
plotted as a function of the product of the atomic numbers A x B, the comparison of 
the charmonium production cross-sections at 450 and at 200 GeV exhibits a remark- 
able feature: the fitted value of a was the same, namely a£o = °- 919 ± °- 015 and 
a 2oo — 0.911 ± 0.034. The remarkable compatibility of the two values justifies a simul- 
taneous fit of the two sets of points, imposing a single a exponent. This global fit leads 
to a J ^ = 0.918 ± 0.015. The observed unexpected agreement between energies and 
reactions is a serious double hint that nuclear dependence is identical for 450 GeV and 
for 200 GeV and that, moreover, charmonium production in light ion-induced reactions, 
namely O-Cu, O-U and S-U, can be considered as normal, namely identical to the one 
observed in p — A interactions. From the global fit can also be determined the ratio be- 
tween the values of <7o for the 200 and 450 data sets. It amounts to 0.38 ± 0.04 and 
results from the changes in both and in the rapidity domain covered by the two sets 
of data. The same factor can be used to rescale the 450 GeV ip' cross-sections to the 
incident momentum and rapidity range of the 200 GeV data. In fact the ratio between 
the if)' and J/tp production cross-sections in p — A collisions seems to be independent 
of y/s, within the energy domain covered by the currently existing measurements [92]. 
The J/ip cross-sections per nucleon-nucleon collision, rescaled if necessary as described 
above, are shown in Fig. 2. The excellent agreement between the 200 and the 450 GeV 
results can be easily judged in the case of p-Cu and p-W collision systems for which both 
measurements exist. 

Within the scenario of a dissociation cross-section, its numerical value can also be ex- 
tracted from the data in the frame of a Glauber calculation. Fig. 3 compares the results, 
rescaled to 200 GeV and to the 0.0-1.0 rapidity interval, with an overall fit which leads 
to a- a bs = 7.1 ± 1.2 mb. The data exhibit a nice overall agreement with the "Glauber" 
calculated values and suggest good compatibility with the model, at their level of preci- 
sion (systematic uncertainties have to be taken into account here). Assuming now the 
validity of the model, the most precise value of <r a b s can be determined using only the 3 
most precise measurements, namely p-C, p-Cu and p-W as obtained at 450 GeV. The 
resulting value obtained from these cross-sections is 5.6 ± 0.4 mb. 

The agreement of the Glauber calculation with the experimental data strongly suggests 
that the A-dependence of J ftp production can be properly accounted for by final state 
nuclear absorption up to and including S-U reactions. In particular, there is no indication 
for any new absorption or dissociation mechanism which would suppress the S-U value 
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Figure 1: The 450 and 200 GeV J/ip cross- 
sections per nucleon (times branching ra- 
tio into dimuons, B^), in the rapid- 
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Figure 2: J/ip cross-sections per nucleon, 
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Figure 3: J/ip cross-sections per nucleon, times b.r., plotted as a function of L. The 
450 GeV values are rescaled to the 200 GeV kinematics. The closed circles are the results 
of the "Glauber" calculation. 



relative to the reference baseline established from the p-A systematics. 

The same experiments mentioned above did provide the corresponding first results for 
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ip' production. Statistical errors, of the order of 10%, were very much too large for 
precise conclusions. In order to minimize systematic effects, the ratios ip' / J/ip measured 
at 450 GeV in p-C, p-Al, p-Cu and p-W can be considered, as taken under identical 
experimental conditions and subject to minimal systematic effects with respect to each 
other. Their values lead to —0.060 ± 0.038 already suggesting, although 

within large errors, a normal absorption for ip' higher than for J / ip . 

In the limited kinematical domain explored by the experiments, the first set of data 
from experiments NA38 and NA51 led, within their uncertainties, to the following strong 
indications : 

• J/ip and ip' are absorbed in nuclear matter. 

• For J/ip, absorption is similar for 450 GeV p-A collisions and for 200 GeV p-A and 
light ion collisions, up to S-U. 

• The absorption is stronger for ip' than for J/ip . 

Much more precise results obtained later by experiment NA50 allowed to confirm these 
results both thanks to significantly increased statistics and to a thorough study of sys- 
tematic uncertainties. 

In order to have full control on systematic effects, Drell-Yan muon pairs, always simul- 
taneously measured with charmonium, are the ideal tool despite their significantly lower 
production cross-section. They allow to use the ratio J/ip / Drell-Yan, indeed with 
increased statistical erors but practically free from systematic uncertainties which are 
usually much more difficult to bring under control. Moreover, they allow with minimal 
theoretical input the study of nucleus-nucleus collisions as a function of centrality. In- 
deed, NA50 has proven experimentally that Drell-Yan production is proportional to its 
theoretically computed value or, equivalently, to the number of nucleon-nucleon collisions, 
from p-p up to Pb-Pb interactions (Fig. 4). It therefore exhibits no nuclear dependence. 
The coherent study of the ratio J/ip / Drell-Yan, both for p-A collisions at 450 GeV (the 
so-called low intensity (LI) although high statistics sample) and for some p-A and S-U 
interactions at 200 GeV, the latter as a function of centrality, leads to another estimate of 
the absorption cross-sections [94]. An independent "Glauber" fit on the S-U ratios gives 
6.3 ± 2.9 mb whereas a simultaneous fit to all the data gives 4.3 ± 0.6 mb, as illustrated 
in Fig. 5. This results further support previous hints that within the errors, no sizable 
additional suppression mechanism is present in S-U collisions with respect to p-A. 

1.3 The First Hints of an Anomaly in Pb-Pb Collisions 

The experimental evidence of an unexpected J/ip suppression resulted from the first 
Pb-Pb data sample, collected in 1995 [95]. It is illustrated on Fig. 6. The anomalous 
character results from the comparison of the J/ip production cross-section per nucleon 
in Pb-Pb interactions with the assumed normal behaviour. The latter, as explained in 
detail above, was established by simple extrapolation of the behaviour obtained from the 
"simultaneous" fit method applied to proton and light ion- induced interactions. 

The plot should be taken as a strong but only qualitative evidence because of the following 
caveat. The experimental reference is based here on the very first p-A and subsequent 
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Figure 4: The Drell-Yan "K-factor", as 
determined experimentally, from p-p up to 
Pb-Pb interactions (from [93]). 
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Figure 6: The J/ip cross-section per nucleon, times b.r., with the 450 GeV results rescaled 
to the same kinematical domain as 200 GeV data, using the "simultaneous" fit method, 
(left) as a function of AxB and (right) as a function of L. The Pb-Pb result, obtained at 
158 Gev, is also rescaled to 200 GeV. 



S-U results, which suffer from both non negligible statistical and also systematic errors. 
The latter have to be taken into account when using results from different setups (NA38, 
NA51). Moreover, it is also based on the implicit assumption contained in the "simultane- 
ous" fit procedure of an energy independent normal absorption cross-section, as suggested 
by the results available at that time. 
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1.4 Anomalous J /if) Suppression in Pb-Pb Collisions 



From the early first indication, it took several years to learn how to make the measurement 
of charmonium suppression in Pb-Pb collisions. The goal was finally reached with three 
different sets of p-A measurements performed from 1996 until year 2000 and two sets of 
Pb-Pb data collected in 1998 and 2000. 

The results of the systematic set of measurements performed at CERN for 450 and 
400 GeV p-A reactions do prove that at these energies, both for J /if) and for if)' : 

• the agreement of the results with the " Glauber model parametrization" is extremely 
good. 

• it becomes excellent when systematic uncertainties are perfectly under control. This 
is the case, for example, when the ratios J/if> / Drell-Yan and if)' / Drell-Yan are 
used instead of the cross-sections themselves, or when absolute cross-sections mea- 
surements become independent from incident beam flux uncertainties (which is the 
case for the NA50 measurements at 400 GeV). 

The two points are illustrated in Figs. 7 and 8 [96,97]. 




Target mass number L (fm) 



Figure 7: The J /if) and if)' cross-sections per nucleon adjusted by the "Glauber model" 
parametrization (left) and compared to the pL parameterization (right) for the p-A data 
collected at 400 GeV in the rapidity range —0.425 < ycM < 0.575. The band (left) repre- 
sents the error associated to the normalization and absorption cross-section uncertainties. 
Within the Glauber model, the fits extrapolation to A=l or to L = (pL) / p =0 leads to 
the value of a . 



1.4.1 J/if> suppression versus the centrality of the collision 

The last sample of Pb-Pb collisions data collected by experiment NA50 could benefit 
from the experience from previous data collections and consequent improvements of the 
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Figure 8: Same as Fig. 7 for the ratios J/ip / DY and ip' / DY. 

detector. It also benefited from overall full coherence between data and simulations, which 
slightly affected normalizations. All the details and results can be found in [98]. 

As displayed in Fig. 9, as a function of the neutral transverse energy Et used here as 
the centrality estimator, the ratio of cross-sections B^aj/^/acY persistently decreases, 
from peripheral to central collisions by a factor of ~ 2.5, showing no saturation in the 
decrease even for the most central collisions. The absorption pattern is here compared 
to the normal nuclear absorption curve as determined from the most recent and accurate 
p-A results obtained in the same experiment at 450 and 400 GeV together with results 
obtained from S-U data collected at 200 GeV. The technique of the "simultaneous" fit 
described in section 3.2 leads to a normal absorption cross-section cr a b s = 4.2±0.4 mb and 
provides part of the rescaling factor needed, from 450/400 GeV to 200 GeV. The advantage 
of using S-U results here is that it leads to minimal uncertainties in the energy rescaling 
factor which is "experimentally" determined, in part at least, as explained before. The 
curve is then further analytically rescaled to 158 GeV under the assumption that cr a b s is 
energy independent. The comparison with the normal absorption curve shows that the 
data behave normally for the most peripheral collisions while increasingly departing from 
this normal behaviour with increasing centrality. 

It has been pointed out that this normal suppression reference could be biased by the use 
of S-U data in its determination since comoving produced hadrons, for example, could 
already affect J/ip production in S-U reactions. A new determination of the reference 
curve has therefore been implemented as detailed hereafter. Use has been made of the 
most precise results on absolute J/ip production cross-sections obtained by NA50 in p-A 
collisions at 450 and 400 GeV together with the results obtained at 200 GeV, again in 
p-A collisions exclusively, by experiments NA38 [99] and NA3 [100]. The Glauber "simul- 
taneous fit" method detailed in section 3.2 leads then to an absorption cross-section 
Cabs = 4.1 ± 0.4 mb which, it is worth noting, is in excellent agreement with the value 
of 4.2 ± 0.4 mb obtained from the ratio of cross-sections J/tp / Drell-Yan. Incidentally, 
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Figure 9: The ratio J/ip / Drell-Yan compared to the normal absorption cross-section 
deduced from p-A and S-U measurements. The inset shows the ratio data / ( normal 
suppression) . 

from the precisely measured ratio of J/ip cross-sections in p-p and p-Pt by experiment 
NA3 at 200 GeV, one gets, from a completely independent determination, using the same 
Glauber approach, an absorption cross-section of o- a b s = 4.1 ± 1.0 mb when restricted 
to the rapidity interval of NA50 Pb-Pb data, which, independently, strongly supports 
the same numerical value for a abs at 450/400 and 200 GeV incident proton energies. As 
detailed in [101], Figure 10 (left and center) illustrates the method and further shows 
that results from O-Cu, O-U and even S-U data exhibit, within errors, a so-called normal 
p-A like behaviour. Figure 10 (right) confirms, that with respect to a pure p-A reference 
curve, J/ip is "anomalously" suppressed in Pb-Pb collisions at 158 GeV. It is worthwhile 
underlining here that the new reference curve for the ratio J/ip / Drell-Yan at 158 GeV 
is partly deduced by analytical rescaling: the J ftp cross-section is rescaled from 200 to 
158 GeV with a Schuler-type formula [97,102] and, on the other side, the appropriate factor 
is applied to the Drell-Yan cross-section to rescale it from 450 to 158 GeV. This pure p-A 
"normal" reference curve has identical shape as the one which was also making use of 
S-U data. It is globally lower by a factor of 0.6% whereas its experimental uncertainty 
is increased by a factor 2. The most recent results obtained from Pb-Pb can now be 
compared with this new exclusively p-A based reference curve. As a function of three 
independently measured quantities tagging the centrality of the collision (see section 2), 
the J/ip survival pattern exhibits very similar trends, as shown in Fig. 11. 

The three independent analyses are displayed together in Fig. 12 after conversion of their 
originally used centrality estimator to the average number of participant nucleons in the 
reaction, N pait . The remarkable agreement between the results of the three analyses, 
underlines the coherence, equivalence and deep understanding of the data provided by 
the three independent detectors used to estimate the centrality of the collision. 

It is worth recalling that E T and the charged particle rapidity density at mid-rapidity 
are directly related to the energy density reached in the collision, through the Bjorken 
formula [109], while the forward energy measured in the ZDC is more strongly correlated 
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Figure 10: J/ip cross-sections per nucleoli in p-A collisions as a function of L. Left: 
Results obtained at 450, 400 and 200 GeV separately. Center: Same results rescaled 
to 200 GeV and compared with results from light nuclei collisions. Right: Same results 
rescaled to 158 GeV and further compared with J /if) production in Pb-Pb collisions. 
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Figure 11: The ratio of cross-sections J ftp / Drell-Yan in Pb-Pb collisions as a function of, 
from left to right, the neutral transverse energy, the "very forward" hadronic energy and 
the charged multiplicity. The curve displays the normal suppression pattern as deduced 
from p-A interactions only. The insets show the ratio data / (normal suppression). 



8 40 



aa 




a(abs) = 4.2 mb (GRV "4 LO) 



50 100 150 200 250 300 350 400 

N__ 




50 100 150 200 250 300 350 400 



Figure 12: The J/ip / Drell-Yan cross-sections ratio as a function of iVp art from three 
analysis originally using the different centrality estimators, compared to (left) and divided 
by (right) the normal nuclear absorption pattern. 
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to the geometry of the collision, and a simple robust estimator of N pait . Note that at 
SPS energies, for a given interacting system, both and are linearly proportional 
to iVp art , as expected in the framework of the wounded nucleon model, up to the most 
central collisions (see, for instance, Refs. [66] and [103]). 

1.4.2 The transverse momentum dependence of J/ip suppression 

The pt spectra of the surviving J / if) 's have been studied by experiment NA50 [104] and, 
in particular, from their most precise sample of Pb-Pb interactions at 158 GeV collected 
in year 2000 [105]. The value of < p\ > for the surviving J/if> 's is plotted in Fig. 13 as a 
function of E^. The overall behaviour exhibits a steady increase with centrality, from the 
most peripheral collisions and up to .Sp— 60 GeV. This value corresponds to L ~ 8.23 fm 
and e = 2.97 GeV/fm 3 . The data also show that, for high centralities, no increase is seen 
up to the most central collisions. 

In order to investigate the underlying mechanism that leads to the observed increase 
of < p\ > for the more peripheral collisions, a comparison is made between different 
colliding systems as shown in Figs. 15 and 16. As a function of L, the same linear 
dependence is observed whatever the colliding system, with a slope apparently indepen- 
dent of the energy of the colliding nuclei. The first observation of this behaviour at 
200 GeV based on results from NA3 [100] and NA38 [106, 107] have been successfully 
interpreted in terms of initial-state parton multiple scattering and analytically described 
by < p\ >=< p\ > pp +a gN x L. 
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Figure 13: The average < p\ > as a func- 
tion of centrality, estimated from E^, for 
the surviving J / ip 's in Pb-Pb collisions at 
158 GeV. 



^-300 
> 





250 



200 



150 



J/ty Pb-Pb 



♦ . ♦ 



♦ + 



+ ♦ 



NA50 2000 (158 GeV/c) 



20 40 60 80 100 120 140 

E T (GeV) 



Figure 14: The transverse mass inverse 
slope parameter T of the surviving J /if} 's 
function of centrality, estimated from 
Et, in in Pb-Pb collisions at 158 GeV. 



All the samples of J /if) considered here are roughly in the same range, namely 
3 < yinb < 4. The values of a g N are nicely compatible for all colliding systems and en- 
ergies. Using the technique of the "simultaneous fit" for the slope and leaving free the 
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Figure 15: The average < p\ > of the sur- 
viving J/ip 's as a function of centrality, 
here L, for p-A (400 GeV), S-U (200 GeV) 
and Pb-Pb (158 GeV) with a "simultane- 
ous" common linear fit. 
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Figure 16: Same as Fig. 15, including 
results on p-A (200 GeV) from NA3 
and NA38, O-Cu and O-U (200GeV) 
from NA38. All results "experimentally" 
rescaled to the kinematical window used 
for Pb-Pb (158 GeV). 



normalization factor to account for the energy change of each set of data leads to the 
following results: 

• < p\ > pp increases linearly with a/s, the total energy in the nucleon-nucleon center 
of mass system. 

• As a function of L, the mean squared transverse momentum of the J/ip exhibits the 
same increasing slope for p-A, O-Cu, O-U, S-U and peripheral Pb-Pb interactions 
as shown in Fig. 16, but saturates for the most central Pb-Pb collisions. 

• The common fit of the increasing slope leads to a g N = 0.081 ± 0.003. 

In an attempt to link this transverse momentum to the rate measurements, it can be 
empirically noticed that, as already discussed above, the measured J/if> suppression rate 
has also the same, so called normal, L dependence from p-p up to S-U reactions. The 
specific anomalous suppression only appears for the most central Pb-Pb reactions, namely 
for L > 7 fm. This peculiar transverse momentum behaviour is also clearly seen from 
Fig. 14 which shows the transverse mass inverse slope parameter T, for the surviving 
J/ip 's, as a function of E^. T is deduced from a fit of the transverse mass distributions 
with the function 1/T • M| • K^Mt/T), where K x is the modified Bessel function. The 
inverse slope, T, is related to the effective temperature of the system in thermal models 
of particle production [107,108]. The measurements show that T increases with centrality 
for peripheral Pb-Pb reactions and then, slightly after anomalous suppression sets in, 
exhibits a trend compatible with no further increase up to the most central collisions. 

The pt dependence of J/ip suppression has been further studied in order to investigate 
what are the pt features of the suppressed J / ip 's, as a function of the centrality of the 
collision. In order to be independent from any external input related to the distribution 
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Figure 17: The ratio of the J /if) p^ distribution in a given bin of centrality to the same 
distribution for the most peripheral bin. With respect to the latter, the low pr part of 
the distributions depletes with increasing centrality. 



of the suppressed J /if) 's, NA50 has limited its study to the ratio of the pt spectra obtained 
in various centrality bins with respect to the spectrum obtained in the most peripheral 
one (a procedure which will be later known as Rcp)- For statistical reasons of this two 
dimensional study, J /if) and Drell-Yan events have been considered within eight different 
centrality bins only. Within each centrality bin, the pt distribution of the J/ if) has been 
studied, normalized to the total number of Drell-Yan events with the same centrality, i.e. 
to the "average number of nucleon-nucleon collisions" in the considered bin. Thus, the 
integral of the pr spectrum in each of the bins is the ratio J /if) / Drell-Yan, as would be 
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plotted, for example, in Fig. 9 dividing the whole centrality range in only 8 bins (instead 
of 11). The ratios Fj and R { are thus defined as: 

dN^/dp T Fi 

F i = Wi and R i = -ET 

iv DY(M>4.2 GeV/c 2 ) ri 

where i is the I th centrality bin, dNj^/dpr is the number of J/ip events of a given 
in centrality bin i and N z DY (m>42 GcV/c 2 ) * s the total number of Drell-Yan events with 
M > 4.2 GeV/c 2 in centrality bin i. 

Figure 17 displays the p T dependence of the seven ratios Ri. It shows that with increasing 
centrality, the lower the the higher the suppression. Moreover and within statistical 
uncertainties, it suggests that, for transverse momenta larger than 2.5 GeV/c, the shape of 
the pt spectrum of J / ip is independent of centrality and very close to the one observed in 
the most peripheral collisions. It should be underlined here that, as all R C p studies, this 
approach does not allow to easily disentangle between transverse momentum spectrum 
modifications due to normal and to abnormal suppression separately. As an illustration, 
a similar study made for S-U reactions which are only subject to normal J/ip suppression 
exhibits nevertheless qualitatively similar features [107]. 

1.4.3 J /if) survival pattern: from p-p to Pb-Pb 

Figure 18 shows the ratio of cross-sections J/ip / Drell-Yan for the collision systems studied 
up to now, from the lightest p-p up to the heaviest Pb-Pb [98]. When necessary, they have 
been reanalyzed in order to ensure a fully coherent data selection and treatment. They 
have also been rescaled to an incident beam momentum of 158 GeV/c when required. 
This rescaling, both for energy and rapidity coverage, has used the "simultaneous fit" 
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Figure 18: The J/ip / Drell-Yan ratio of cross-sections vs. L, for several collision systems, 
compared to (left) and divided by (right) the normal nuclear absorption pattern. All data 
are rescaled to 158 GeV/nucleon. 
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Figure 19: Same as Fig. 18 for S-U and Pb-Pb collisions only, as a function of energy 
density. The absorption curves for S-U and Pb-Pb in the left panel of the figure are 
slightly different because the relation between energy density and L (obtained from a 
Glauber calculation) depends on the colliding nuclei. 



method to bring the results obtained at higher momenta down to 200 GeV/c and an ana- 
lytical calculation to further bring the 200 GeV results down to 158 GeV. Figure 19 shows 
the same ratio of cross-sections J ftp / Drell-Yan obtained in S-U and in Pb-Pb collisions 
as a function of e, the energy density averaged over the whole transverse area of the 
collision. The latter is obtained with the Bjorken [109] formula from the total transverse 
energy deduced from the measurement of its neutral component in the electromagnetic 
calorimeter [113]. It shows that the departure from the normal nuclear absorption curve 
sets in for energy densities around 2.5 GeV/fm 3 , just above the values reached in the most 
central S-U collisions. 

The numerical values of the different variables used as centrality estimators are given in 
Table 1 for the Et bins chosen above for Pb-Pb collisions at 158 GeV. They have been 
calculated with the Glauber model and assuming an initial formation time r = 1 fm/c. 
The values for S-U collisions at 200 GeV are given in Table 2. They include, in the 
errors, the resolution of the electromagnetic calorimeter. The authors [113] consider that 
the relative scale between S-U and Pb-Pb is robust, while the absolute values could be 
affected by a 20%-25% systematic uncertainty. 



2 Features of ip' Suppression at SPS Energies 

As known since long, ip' production is also affected by nuclear effects in p-A collisions. 
This is now quantitatively supported with precise results from NA50 based on two differ- 
ent samples of data collected at 450 and 400 GeV with incident protons on Be, Al, Cu, Ag, 
W and Pb targets [96,97,110,111]. The results, obtained with a Glauber parametrization 
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Bin 




(GeV) 


-Vpart 




b (fm 


) 


L (fm) 


e (GeV/fm 3 ) 




range 


average 


average 


rms 


average 


rms 


average 


rms 


average 


rms 


1 


3-15 


10.6 


34 


13 


11.8 


0.7 


4.44 


0.72 


1.15 


0.34 


2 


15-25 


20.4 


70 


13 


10.3 


0.5 


5.94 


0.42 


1.86 


0.20 


3 


25-35 


30.3 


104 


14 


9.2 


0.4 


6.84 


0.33 


2.27 


0.15 


4 


35-45 


40.2 


138 


16 


8.2 


0.4 


7.51 


0.27 


2.60 


0.13 


5 


45-55 


50.2 


172 


17 


7.3 


0.4 


8.02 


0.23 


2.87 


0.11 


6 


55-65 


60.1 


206 


18 


6.5 


0.4 


8.43 


0.20 


3.06 


0.09 


7 


65-75 


70.1 


240 


19 


5.6 


0.5 


8.76 


0.17 


3.27 


0.08 


8 


75-85 


80.1 


274 


20 


4.8 


0.5 


9.02 


0.15 


3.46 


0.07 


9 


85-95 


90.1 


308 


21 


3.9 


0.6 


9.22 


0.12 


3.66 


0.06 


10 


95-105 


100.0 


341 


20 


2.8 


0.7 


9.38 


0.11 


3.83 


0.05 


11 


105-150 111.5 


367 


15 


1.7 


0.7 


9.48 


0.06 


3.89 


0.04 



Table 1: Centrality classes for Pb-Pb interactions at 158 GeV based on the transverse 
energy measurement. For each class are listed the E T range and weighted average, together 
with the average and rms values of N part , b, L and e. 



Bin 




(GeV) 


-Vpart 




L (fm) 


e (GeV/fm 3 ) 




range 


average 


average 


rms 


average 


rms 


average 


rms 


1 


13-28 


22.1 


31.3 


9.1 


4.37 


0.54 


1.04 


0.22 


2 


28-40 


34.3 


50.0 


9.7 


5.35 


0.43 


1.46 


0.18 


3 


40-52 


46.3 


66.7 


10.8 


6.04 


0.41 


1.76 


0.16 


4 


52-64 


58.3 


83.5 


11.7 


6.65 


0.41 


2.01 


0.15 


5 


64-76 


70.2 


98.8 


11.0 


7.17 


0.38 


2.22 


0.13 


6 


76-88 


81.7 


109.0 


8.2 


7.53 


0.29 


2.38 


0.09 



Table 2: Same as previous table, for S-U collisions at 200 GeV. 



separately at 450 and 400 GeV, are quite compatible. Their combination leads to an 
absorption cross-section of 8.3±0.9 mb from the values of production cross-sections per 
nucleon and 7.7±0.9 mb from the ratios ip' / Drell-Yan, with minimal systematic uncer- 
tainties in this case [96]. These values can be compared to 4.18±0.35 mb as obtained for 

The ip' study [112] made with Pb-Pb interactions at 158 GeV is based on only 1285 
■ip' events. It shows that both with respect to Drell-Yan, and also with respect to the 
J/ip , ip' production decreases with increasing centrality, as displayed in Figs. 20 and 21. 
It is worthwhile noting here that the Drell-Yan mass range chosen here to normalize the 
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Figure 20: The ratio of cross-sections 
B' cr(^')/cr(DY) as a function of E T , with 
DY in the mass range 4.2-7.0 GeV/c 2 . 
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Figure 21: The ratio B'^a(ip')/ fi w a(J/^) 
as a function of Et- 



ip' yield is 4.2-7.0 GeV/c 2 . Scaling up the results to the mass range used in the J /if} sup- 
pression studies, namely 2.9-4.5 GeV/c 2 , requires to apply the factor 7.96, which is the 
ratio between the Drell-Yan cross-sections in the two mass domains. In order to investi- 
gate whether ip' production exhibits any abnormal suppression pattern in nucleus-nucleus 
collisions, the production cross-sections per nucleon-nucleon collision (or, equivalently, 
normalized to Drell-Yan ) are compared, after appropriate rescaling to the same kinemat- 
ical domain, between p-A, SU and Pb-Pb collisions [96,97]. 



> 

D 



e 

i: 



All data rescaled to 158 GeV 



A NA50 p-A 450 GeV, HI 96/98 
* NA50 p-A 450 GeV , LI 98/00 
o NA50 p-A 400 GeV, HI 2000 



10 




£1.8 



--1.6 

"D 
CD 

^1.4 



E1.2 
1 

0.8 
0.6 
0.4 
0.2 



1r ^ " r 



i Pb-Pb 

■ S-U 

a p-A HI 96/98 

* p-A LI 98/00 

o p-A HI 2000 



tt, 



9 10 
L(fm) 



Figure 22: The ip'/ Drell-Yan ratio of cross-sections vs. L, for several collision systems, 
compared to (left) and divided by (right) the normal nuclear absorption pattern. All data 
are rescaled to 158 GeV/nucleon. 



The left panel of Fig. 22 shows the ratio B' a{ip') / cr(DY) as a function of L. The measured 



47 



B 
o 

©1.8 

Q- 

X 

^1.6 

"8 

=1-4 

w 
nj 

Jj1.2 
1 

0.8 
0.6 
0.4 
0.2 




■ ■ Pb(208x 158GeV/c)-Pb 
-• S(32 x 200 GeV/c)-U 

■ * p(450 GeV/c, Ll)-A 
- A p(450 GeV/c, Hl)-A 

■ T p(400 GeV/c, VHI)-A 

3 12 3 4 



□ Pb(208x 158G8V/o)-Pb 

O S(32 x 200 GbV/c)-U 

* p(450 GeV/c, Ll)-A 

A p(450 GeV/c, Hl)-A 

« p(400 GeV/c, VHI)-A 



♦tt, 



9 10 
L(fm) 



Figure 23: Same as Fig. 22 (right) for 
J/ip and tp' on same L scale. 
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Figure 24: The ratio B'^a(^')/a (DY) as 
function of the energy density. 



suppression patterns suggest the following features: 

• A fair agreement with exponential behaviours, with two different regimes, one for 
proton and a different one for ion-induced reactions. 

• A similar centrality dependence for S-U and Pb-Pb interactions. 

Using the approximate exponential parametrization to describe ip' absorption as a function 
of L, the fit of the data leads to an absorption cross-section of 7.3 ± 1.6 mb in p-A 
collisions, while a much higher value, 19.2 ± 2.4 mb, is obtained for ion-ion collisions (S-U 
and Pb-Pb fitted simultaneously). The right panel of Fig. 22 shows, as a function of L, the 
ratios between the measured ip' yields, normalized to Drell-Yan, and the corresponding 
expected "normal nuclear absorptions", as extrapolated from p-A measurements. The 
difference between the suppression patterns for J/ip and tp' is shown in Fig. 23. The 
ip' suppression for S-U and Pb-Pb collisions is shown in Fig. 24 as a function of the 
energy density which is calculated as described in detail in [113]. The numerical values 
of L, iVp art and the energy density e used for the ip' studies in Pb-Pb collisions are given 
in Table 3. 



3 More Results from SPS and RHIC 

After the systematic studies carried out by experiments NA38, NA50 and NA51, two 
other experiments were aimed at extending the knowledge on charmonium suppression in 
ultrarelativistic heavy ion interactions. 

At the CERN-SPS, in a fixed target experiment, the NA60 collaboration uses basically 
the same muon spectrometer as NA50 [88] except for the centrality estimate which is 
based exclusively on the very forward calorimeter measurements. The NA50 detector is 
upgraded with a silicon pixel telescope, located in the target region, which leads to an 
extremely precise determination of the origin of the measured tracks and, consequently, 
to a better dimuon mass resolution. Data have been collected to study In-In collisions 
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Bin 




(GeV) 


^Vpart 




L (fm) 


e (GeV/fm 3 ) 




range 


average 


average 


rms 


average 


rms 


average 


rms 


1 


3-20 


13.9 


44.6 


17.7 


4.90 


0.84 


1.24 


0.37 


2 


20-35 


28.2 


96.7 


18.0 


6.65 


0.44 


2.04 


0.20 


3 


35-50 


43.0 


147.1 


19.3 


7.65 


0.31 


2.53 


0.14 


4 


50-65 


57.8 


197.7 


20.6 


8.34 


0.24 


2.89 


0.11 


5 


65-80 


72.7 


248.6 


22.0 


8.83 


0.18 


3.19 


0.09 


6 


80-95 


87.6 


299.2 


23.0 


9.17 


0.14 


3.52 


0.07 


7 


95-150 


105.9 


353.3 


22.4 


9.43 


0.10 


3.76 


0.06 



Table 3: Mean free path crossed by the cc pair inside the nucleus, number of participating 
nucleons and energy densities for each Pb-Pb E T range. 

with an incident In beam of 158 GeV/c [114] and also p-A reactions with incident proton 
beams of both 400 and 158 GeV/c. 

At the RHIC collider at BNL, the PHENIX experiment studies p-p [115,116], d-Au [117, 
118], Cu-Cu [119] and Au-Au [120,121] collisions at much higher energies, namely 200 GeV 
in the c.m.s. and, consequently, at higher energy densities than those reached at the 
CERN-SPS. It makes use of a detector specifically designed to study heavy ion collisions 
at RHIC [122]. 

Some of the results obtained by these experiments in their studies of J / ip suppression are 
summarized hereafter. 

3.1 J/ip Suppression in In-In Collisions at 158 GeV 

In a first analysis strategy, experiment NA60 follows the approach adopted by experi- 
ments NA38/NA50 and makes use of the ratio J ftp / Drell-Yan. This ratio is free from 
systematic uncertainties and provides its own absolute normalization, without any ex- 
ternal input. It thus leads to the most robust estimate of the J/ip yield normalized to 
the number of nucleon-nucleon collisions, through the experimentally measured Drell-Yan 
events. Unfortunately the number of Drell-Yan events in NA60 is extremely small and 
allows, in a centrality dependent study, for only 3 bins in centrality. In order to identify 
any anomalous behavior, the results are compared with normal nuclear absorption as de- 
termined experimentally by experiments NA38/NA50. The results can then be expressed 
in terms of the ratio "measured/expected", where "expected" refers to normal nuclear ab- 
sorption corresponding to an absorption cross-section of 4.18±0.35 mb. The 3 measured 
ratios are plotted in Fig. 25 as a function of A^ part . 

In a second approach, the measured " unnormalized" J/ip yield is directly compared to 
the analytically calculated centrality distribution of J/ip . The latter is derived in the 
frame of the Glauber model using the same absorption cross-section of 4.18 mb. Relative 
normalization between the data and the reference curve is obtained by requiring that the 
ratio of the centrality integrated distributions is equal to the centrality integrated ratio 
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J/if> I Drell-Yan. Fig. 26 shows both the measured and calculated distributions as a 
function of centrality, estimated here from the forward energy of the spectators nuclei. 
Their ratio is plotted in Fig. 25 as a function of A^ part . 

The method does not suffer any more from the limited number of Drell-Yan events, makes 
use of the statistical power of the J / if> sample and reaches unprecedented small statistical 
uncertainties. It should be noted here, nevertheless, that the same method could have 
been applied by experiment NA50 which, indeed and as compared to NA60, had more than 
a factor 3 larger J /if) sample of events. The NA50 experiment made the choice of studying 
the J /ip suppression pattern in Pb-Pb collisions through the robust, unbiased J /if) / Drell- 
Yan ratio which is free, at first order, from all potential systematic experimental effects 
which cancel out in the ratio. It is also free from analytical model dependent descriptions 
of the detectors. It therefore minimizes systematic uncertainties, sometimes difficult to 
find and evaluate. NA50 discarded this type of analysis from its final results which are 
thus dominated, by far, by pure statistical uncertainties in the Drell-Yan sample of events. 
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Figure 26: The J /if) -Ezdc distribu- 
tion (circles), compared with expectations 
from nuclear absorption (line). 



Figure 25: Centrality dependence of the 
J /if} suppression in In-In collisions for the 
ratio oj/^/ctdy (star symbols) and for the 
J /if} absolute yield (circle symbols). Also 
shown are the common global systematic er- 
ror (left box) and the relative point to point 
uncertainties. 

The results of experiment NA60 can be compared to those of NA50 shown in Fig. 12 
which superposes the results obtained, separately, through the three independent cen- 
trality estimators. The suppression patterns in Pb-Pb and In-In are qualitatively similar 
and compatible, thanks to the errors quoted by NA50. In some more detail, nevertheless, 
some differences can be noticed in this comparison. 

The results on In-In do not exhibit a solid evidence of normal behaviour for peripheral 
collisions. There is one single centrality point with a ratio "measured/expected" amount- 



50 



ing to 0.99±0.06, which is somewhat insufficient to support a constant trend close to 
unity in the vicinity of N paxt ~ 50. On the other hand, the Pb-Pb results do exhibit 
a normal behaviour extending over 2 or 3 centrality bins and up to about 76 iV part . 
Taking into account the resolution of the centrality measurement, In- In collisions would 
depart from normal behaviour for iV part ~ 80 iV part corresponding to an energy density 
of 1.5 GeV/fm 3 . For Pb-Pb collisions, the departure is detected at iV part ~ 100 iVp art 
when using the same centrality detector and at iVp art ~ 120 iVp art when using the most 
precise centrality estimator, which corresponds to 2.25 and 2.45 GeV/fm 3 respectively. 

The In-In and Pb-Pb suppression patterns themselves show different trends. Pb-Pb ex- 
hibits a steady decreasing trend, in particular for A^ part > 150, as opposed to a kind of 
an unexpected rise, starting for A^ part > 140, in the case of In-In. This rise could perhaps 
accommodate, within errors, some flat behaviour. 



3.2 J/ip Suppression in A-A Collisions at y/s = 200 GeV 

The PHENIX experiment at RHIC measures J /if) production at y^s = 200 GeV, for 
rapidities \y\ < [1.2,2.2] and also \y\ < 0.35. The study is made in terms of the "nuclear 
modification factor", -Rab, defined as: 

(da/dy) A - B 



Ncou (da/dy) p _ p 

when referring, for example, to the reaction of nucleus A with nucleus B. 

Fig. 27 shows the ratio -Raa for J/ip 's of low rapidity in the collision c.m.s. obtained by 
the PHENIX experiment for Au-Au and Cu-Cu collisions. The same plot shows the results 
obtained with other colliding systems at lower energies, derived from the measurements 
performed at CERN. The ratio -Raa (or -Rab) is plotted as a function of the iV par t value 
corresponding to the selected centrality bins made in the different analyses. It exhibits, 
within errors, an in-discriminable pattern for all the colliding systems, including S-U and 
Pb-Pb. On the other hand and from different type of studies, as detailed above, it is 
known that S-U behaves normally as well as the 2 more peripheral points of Pb-Pb which 
are plotted in the figure, whereas abnormality sets in for iV par t > 100 in the case of 
Pb-Pb, at least. It therefore appears that this kind of representation, namely -Raa as a 
function of iV part , is not the most sensitive in order to experimentally detect, without 
external theoretical inputs, any abnormal J/ip suppression. 

Fig. 28 shows that the nuclear modification factor is a function of the J/x/j rapidity for 
Au-Au collisions. Indeed, this feature was already visible in fixed target p-A experiments 
with 200 and 800 GeV/c incident protons by experiments NA3 [100] and E866 [124]. 
The understanding of such an observation on Au-Au at 200 GeV requires again complete 
experimental knowledge of the reference for the two different kinematic windows, under 
conditions excluding the formation of a new state of matter (see section "Discussion" 
below) . 

The J/ip suppression has also been studied by PHENIX as a function of the transverse 
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Figure 27: J/ip nuclear modification factors 
for Au-Au, Pb-Pb, In-In and S-U colliding 
systems at their respective energies (200, 19 
and 17.3 GeV ) function of the number 
of participants, N pSuTt . 
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Figure 28: (a) The J ftp nuclear modifica- 
tion factor for Au-Au c\S cX function of N, 
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for central and forward rapidities, (b) Ra- 
tio of forward/mid rapidity J/ip Raa vs. 
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Figure 29: The surviving J/ip average 
< p\ > as a function of iVp art for Au-Au, 
as measured by PHENIX at mid (left) and 
forward (right) rapidity. 
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Figure 30: Same as fig. 29, as a function 
of L. 



momentum. Fig. 29 shows, as a function of A^ part , the values of < p\ > for various collid- 
ing sytems [119]. Fig. 30 [123] uses instead the variable L in order to allow a comparison 
with the results obtained at the SPS, successfully parametrized as a function of the path 
traversed by the J/ip in nuclear matter (see Figs. 15 and 16). The PHENIX measure- 
ments are affected by very large uncertainties. Within their limited accuracy, they also 
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suggest a linear behaviour as a function of L, the same for all the systems which are taken 
here at same colliding energy. The numerical values of the corresponding fitted slopes, 
namely 0.012 ± 0.046 and 0.093 ± 0.035 for central and forward rapidity respectively can 
be compared with the accurate value measured at the SPS by the NA38/NA51/NA50 ex- 
periments which amounts to 0.081 ±0.003 in significantly different kinematical conditions. 
Much more accurate measurements are obviously needed at RHIC in order to become a 
quantitative probe for the simplistic interpretation of the lower energy results. 

4 Discussion and Evaluation 

Future experimental progress could consolidate what appears already today like well es- 
tablished experimental evidences. The results detailed above do raise, nevertheless, two 
important points which are discussed hereafter. 

• Anomalous J/ip suppression in Pb-Pb collisions at ^/s = 17.2 GeV. 

The anomalous character of the observed J/ip suppression is the first point to be 
clarified if charmonium suppression is the signature of some new physics specific 
to high energy nucleus-nucleus collisions. For the moment, an evidence for this 
abnormal behaviour has been deduced from the systematic studies carried at the 
SPS by experiments NA38/51/50. In order to define the reference curve attached 
to " normality" , the experiments have made use of a set of systematic measurements 
made on p-induced reactions at 450, 400 and 200 GeV/c. It has been found that 
also J/if) 's produced in O-Cu, O-U and S-U collisions at 200 GeV behave normally, 
i.e. like in p-A reactions, for the same rapidity range in the lab (or rather target) 
reference system. It is with respect to this experimental reference, appropriately 
rescaled to 158 GeV under the explicit assumption that nuclear absorption of J ftp is 
beam-energy independent from 450 down to 158 GeV, that the abnormal character 
of J/ip suppression has been established for mid-central and central Pb-Pb reactions. 
The validity of this crucial assumption is supported by two independent fixed target 
experiments. 

1. Experiment E866 

At Fermilab, experiment E866 [124] has provided cross-section ratios measure- 
ments on p-Be and p-W at 800 GeV/c (a/s = 39 GeV). The measured ratio 
of cross-sections remains constant within the range —0.10 < Xf < 0.25 or, 
equivalently 3.13 < y\ ah < 4.95. When parametrized with the Glauber model, 
it leads to an absorption cross-section of 2.83 ± 0.77. This numerical value 
is therefore rapidity independent within the explored y^ range which largely 
overlaps with the range 3.0 < t/i a b < 4.0 explored by the NA38/NA51/NA50 
experiments. The latter find 4.18±0.35 mb, show that nuclear absorption in- 
deed accounts for the measurements and, furthermore, that the absorption 
cross-section is rapidity independent within the covered range [110,111]. 

2. Experiment NA3 

The most precise result on the topic is found in CERN experiment NA3 [100, 
125]. Their data lead to a precise value of the ratio of J/ip production cross- 
sections measured directly in p-p and in p-Pt collisions in a single simultaneous 



53 



beam exposure at 200 GeV or y/s = 19.4 GeV. The Glauber model descrip- 
tion allows to derive the nuclear absorption cross-section, with appropriate 
treatment of errors, in the yi a b rapidity range where NA3 shows it stays con- 
stant, namely 3.0 < y^ < 4.0. The resulting numerical value is 4.1±1.0 mb. 
This value is surprisingly close to the NA50 determination at 400/450 GeV, 
namely 4.18±0.35 mb. The results obtained for p-a reactions at 800, 450/400 
and 200 GeV/c suggest that the J/ip absorption cross-section has no or little 
dependence on y/s within this -Ebeam range, when considered in the same y^ 
rapidity interval in which it stays constant. They thus support the assumption 
made by NA50 when rescaling the normal absorption reference curve down to 
158 GeV. 
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Figure 31: -Raa for Au-Au collisions compared to a band of theoretical curves for the 
f^abs values found to be consistent with the d-Au data. EKS and NDSG shadowing are 
included for midrapidity (left). Same for forward rapidity (right). 



• Is there anomalous J / ip suppression in Au-Au collisions at RHIC ? 

In order to be able to detect an abnormal J/ip suppression from an experimental 
point of view and with minimal external input, appropriate experimental references 
are mandatory, measured under conditions such that no new physics can reasonably 
be expected. The PHENIX experiment has made use of the presently available d-Au 
data to establish the normal suppression or "cold nuclear matter effects" pattern ex- 
pected for the Au-Au measurements at RHIC [126]. Figs. 31 and 32 show the values 
of Raa as a function of centrality as measured in Au-Au collisions compared to the 
pattern extracted from d-Au measurements, taking also into account the shadowing 
effects expected at RHIC. PHENIX conclusion is that "Neither the predictions of 
cold nuclear effects in heavy ion collisions based on fitting of the d-Au data with 
theoretical curves nor those obtained directly from the d-Au data points are well 
enough constrained to permit quantitative conclusions about additional hot nuclear 
matter effects" . Much more precise reference data on d-Au collisions are needed to 
reach a solid conclusion. 
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Figure 32: -Raa for Au-Au collisions compared to predictions of the data driven method 
constrained by the Rd-Au as a function of collision centrality, for midrapidity (left) and 
forward rapidity (right). 



5 Summary of the Experimental Status 

At the eve of the start of a new generation of experiments aimed at further studies of 
quarkonium abnormal suppression as a signature of QGP formation, the present experi- 
mental situation is summarized hereafter. 

1. The production of J /if) suffers from cold matter nuclear effects also called normal nu- 
clear absorption. This absorption can be interpreted as resulting from interactions 
of a given cc state while propagating in normal cold nuclear matter. Quantitatively, 
the corresponding absorption cross-section has been precisely measured in p-A in- 
teractions up to p-Pb, at y/s = 29.0 and 27.4 GeV. For a J/ip with rapidity in the 
range 3.0 < yi a b < 4.0. with respect to the traversed nuclear matter, the best de- 
termination of the absorption cross-section leads to cr a b s J ^ = 4.2 ± 0.35 mb. The 
values derived from data collected at = 39 and 19.3 GeV in the same rapidity 
range are in in very good agreement with this devoted measurement. 
Measurements performed on light nuclei interactions, namely O-S, O-U and S-U 
at 19.3 GeV show excellent compatibility with the above quoted cross-section and, 
therefore, J /if) production in such reactions only suffers from normal nuclear effects. 
There is thus compelling evidence that the value of cr a b s in the same rapidity window 
stays the same for -Js = 17.3 

2. The production of if>' also suffers from cold matter nuclear effects. Systematic mea- 
surements performed on p-A reactions at 29.0 and 27.4 GeV lead to a cross-section 
Cabs ^ = 7.7 ± 0.9 mb, significantly higher than a a bs J ^ ■ 

3. From comparison with appropriate references deduced from the measurements on p 
and light ion-induced collisions, J/if> production in Pb-Pb interactions at 17.3 GeV 
exhibits only normal nuclear effects up to an energy density of e = 2.5 GeV/fm 3 . 
An abnormal absorption appears above this value of e and significantly increases 
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with increasing energy density, up to an abnormal suppression factor of 0.6, for the 
most central Pb-Pb collisions. 

4. In contrast with the J /ip patterns, ip' production exhibits an abnormal suppression 
starting with S-U reactions at 19 GeV. The observed suppression pattern is in very 
good agreement with the one exhibited by Pb-Pb interactions at 17.3 GeV. They 
both lead to an absorption cross-section of <7 a b s ^ = 19.2 ± 2.4mb as compared to 
7.3 ± 1.6 for normal ip f absorption. Abnormal suppression for ip' starts above an 
energy density e = 1 GeV/fm 3 . 

5. Data obtained from In-In collisions at 17.3 Gev tend to qualitatively support the 
suppression pattern measured in Pb-Pb collisions. 

6. For the moment being, there is no clear abnormality in Au-Au collisions at 200 GeV. 
Further measurements to precisely quantify cold nuclear effects are a must. 

7. J/tfj absorption, as observed in Pb-Pb collisions, affects only low p T J /if) 's for 
peripheral collisions. Higher J/ ip 's are more and more suppressed with increasing 
centrality of the collision. This feature requires a robust experimental reference in 
order to allow model-independent conclusions. 

8. The surviving J/ip 's in Pb-Pb collisions exhibit an average < p\ >, or equivalently, 
a temperature T e s which starts by increasing with increasing collision centrality. It 
saturates for energy densities above ~ 2.8 GeV/fm 3 and up to 3.9 GeV/fm 3 , the 
average value of the most central bin. 



Outlook 

The theoretical analysis of the in-medium behaviour of quarkonia has greatly advanced 
in the past decade. Potential model studies based on lattice results for the colour-singlet 
free energy as well as direct lattice calculations appear to be converging, and within a 
few more years, the dissociation temperatures for the different quarkonium states will be 
calculated precisely. Through corresponding calculations of the QCD equation of state, 
these temperatures provide the energy density values at which the dissociation occurs. In 
statistical QCD, quarkonia thus allow a spectral analysis of the quark-gluon plasma. 

The application of this analysis in high energy nuclear collisions is so far less conclusive. 
There exists a wealth of data from different interactions at the SPS, and first RHIC results 
are now also available. These results indicate the production of a hot, dense, strongly 
interacting medium. It is not yet clear, however, to what extent this medium is indeed 
the quark-gluon plasma studied in statistical QCD. Is it possible to identify dissociation 
onsets for different charmonium states? Does the initial (non-thermal) overabundance of 
charm survive a subsequent thermalization? Is the apparent J/ip suppression saturation 
from SPS up to fairly central RHIC collisions the survival of the directly produced J/^'s as 
seen in lattice QCD studies, or is it the production of additional J/^'s through secondary 
cc pairings? These questions can only be answered by further experiments, ideally in all 
the possible energy ranges. 
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Future nucleus-nucleus experiments at the CERN-LHC will extend the field to much 
higher energy densities; but the precise experimental measurement of cold matter nuclear 
effects in this new energy domain will still remain the prime condition for the understand- 
ing and interpretion of the results. 
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